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The influence of randomly distributed point impurities and planar defects on order and transport 
in type-II superconductors and related systems is studied. It is shown that the Bragg glass phase 
is unstable with respect to planar defects. Even a single weak defect plane oriented parallel to the 
magnetic field as well as to one of the main axis of the Abrikosov flux line lattice is a relevant 
perturbation in the Bragg glass. A defect that is aligned with the magnetic field restores the flux 
density oscillations which decay algebraically with the distance from the defect. The theory exhibits 
striking similarities to the physics of a Luttinger liquid with a frozen impurity. The exponent for the 
flux line creep in the direction perpendicular to a relevant defect is derived. We flnd that the flux 
line lattice exhibits in the presence of many randomly distributed parallel planar defects aligned to 
the magnetic field a new glassy phase which we call planar glass. The planar glass is characterized by 
diverging shear and tilt moduli, a transverse Meissner effect, resistance against shear deformations. 
We also obtain sample to sample fluctuations of the longitudinal magnetic susceptibility and an 
exponential decay of translational long range order in the direction perpendicular to the defects. The 
flux creep perpendicular to the defects leads to a nonlinear resistivity p{J ^ 0) ~ exp[— ( Jd/ J)'*/^]. 
Strong planar defects enforce arrays of dislocations that are located at the defects with a Burgers 
vector parallel to the defects in order to relax shear strain. 

PACS numbers: 74.25.Qt, 71.55.Jv, 74.62. Dh,64.70.Rh 



I. INTRODUCTION 

Type-I superconductors are both perfect conductors 
and perfect diamagnets. In type-II superconductors the 
perfect diamagnetism is reduced, an external magnetic 
field penetrates the sample above the lower critical field 
Hci in the form of magnetic flux fines (FLs)i. A trans- 
port current will then lead to a motion of the FLs, yield- 
ing a linear resistivity p « pnB / Hc2 in disorder-free sam- 
ples. Here B denotes the magnetic induction and Hc2 is 
the upper critical field^. At i? = Hc2 diamagnetism dis- 
appears completely and p reaches the resistivity p„ of the 
normal state. In order to recover the desired property of 
a dissipation-free flow, FLs have to be pinned. Point de- 
fects, such as vacancies or interstitials, are one type of 
pinning source. In high-Tc materials point impurities are 
almost always present due to a non-stoichiometric com- 
position of most materials. Impurity pinning leads to a 
zero linear resistivitj*^. However, thermal fluctuations al- 
low for FL creep, resulting in a non-zero nonlinear resis- 
tivity of the form p{J) ~ exp[— ( Jp/ J)^] where the creep 
exponent is /i = 1/2 for point impurities^. J( ^ Jp) 
denotes the current density and Jp depends on B, tem- 
perature T, concentration and strength of the pinning 
centers as well as on properties of the material through 
the superconductor coherence length ^ and the penetra- 
tion length A. This FL creep law is closely related to 
the order of the flux line lattice (FLL) in the presence of 
point pinning centers. 

The order of the FLL was a puzzle for a long time. 
Larkin concluded in 1970 that randomly distributed im- 
purities destroy the long range order of the Abrikosov 



lattice^. Only much later it was realized that the effect 
of impurities is weaker, resulting in a power law decay 
of translational order of the FLs in the so-called "Bragg 
glass" phaseii^iii^iaiia. 

More effective pinning sources can further suppress the 
nonlinear resistivity. One example are columnar defects, 
produced by heavy ion radiation, that have been consid- 
ered by Nelson and Vinokur— li^. These authors mapped 
the physics of FLs onto the problem of the localization 
of bosons in two dimensions where FLs play the role of 
world lines of the bosons^. At low temperatures they 
found strongly localized FLs at the columnar defects, 
forming a "Bose glass" phase. Thermally activated hop- 
ping of noninteracting FLs in the limit J — > leads to the 
creep exponent p. = 1/3, while FL interactions yield the 
increased creep exponent p = 1—. The transport in this 
regime closely resembles the variable range hopping of 
electrons in two dimensional disordered semiconductors. 
This picture is expected to be valid for weak enough ap- 
plied magnetic fields, such that the density of defects is 
bigger than the FL density. For a larger magnetic field, 
Radzihovsky— argued that the Bose glass coexists with 
a resistive liquid of interstitial FLs which upon cooling 
freezes into a weakly pinned Bose glass. For asymptoti- 
cally weak currents, the creep of FL bundles determines 
the nonlinear resistivity and p = lis the creep exponent^. 

In this work we consider planar defects like twin 
boundaries from which even stronger pinning can be ex- 
pected. Twin boundaries are ubiquitous in supercon- 
ducting YBa2Cu307_a; and La2Cu04 where they are 
needed to accommodate strains arising from tetragonal 
to orthorhombic transition as a result of oxygen vacancy 
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FIG. 1: Schematic phase diagram of disordered flux line lat- 
tices resulting from impurities, columnar and planar defects 
of concentration riimp, tied and rip^, respectively. The stabil- 
ity of the phases with respect to different kinds of disorder is 
indicated by arrows, fx denotes the creep exponent. 



ordering and due to rotation of the CuOg octahedra, 
respectiveljii^. Twin boundaries occur frequently with 
the same orientatio n^^i^^i^^ or in orthogonal families of 
lamellas (" colonies" They can be regularly dis- 

tributed with rather fixed spacing or with large variations 
in the spacingi^. The mean distance of the defect 
planes varies between 10 77,77 ^15, le^is ^^^^ l^rnM^. 

The common feature of all of the above mentioned de- 
fects is that they lead to FL pinning, but what distin- 
guishes them is the nature of the pinned phase. In con- 
trast to point disorder, which promotes FL wandering, 
planar defects inhibit wandering and promote localiza- 
tion. Pinning of individual FLs by columnar as well as 
planar defects in the presence of bulk point disorder has 
been investigated in the past^i^^i^i^l. The competition 
between a planar defect and point impurities in three- 
dimensional systems, for a single FL, leads to localiza- 
tion of the FL at all temperatures^ii^iS^. The influence 
of many parallel defect planes on the creep of a single 
FL perpendicular to the planes has been studied at low 
temperatures when the FL spacing exceeds the average 
spacing between the planesi^. 

The main focus of this paper is correlated disorder 
in the form of planar defects. Some of the results of 
this paper have been published earlie r^^i^^ . Here we give 
additional results and present more detailed derivations. 
First, we discuss the influence of a single planar defect on 
the stability of the Bragg glass phase. Then we explore 
the effect of many defect planes on the FLL. We find that 
the necessary condition for a planar defect to become a 
relevant perturbation is that it is oriented parallel to the 
magnetic field. In this case, its influence on the Bragg 
glass phase can be characterized by the value of a sin- 
gle parameter g = ^ri{a/£)'^ which depends both on the 
exponent t] describing the decay of the positional corre- 
lations in the Bragg glass phase and on the orientation of 
the defect with respect to the FLL. a and i are the mean 
spacing of the FLs in the absence of the defect plane and 
the distance between lattice planes, of the Abrikosov lat- 
tice, parallel to the defect, respectively. A weak defect 
turns out to be relevant if g < 1, i.e., if it is parallel to 
one of the main crystallographic planes of the FLL. 

The FL density averaged over point impurities shows 



periodic order with an amplitude decaying as a power 
law with the normal distance to the defect plane. For 
a relevant defect on scales larger than L^,, the expo- 
nent controlling the power law is g. For the definition 
of Lu, see Eq. (p8)) below. For g > 1 a. weak defect 
is irrelevant (in the sense of rcnormalization group) and 
the density profile decays faster, with a larger exponent 
2(7— 1 > 1, on scales larger than the positional correlation 
length La- For a weak defect tilted against the applied 
magnetic field we find that FL density oscillations de- 
cay exponentially fast. We investigate also the dynamics 
of FL bundles perpendicular to the relevant defect for 
small current densities J — > 0. The nonlinear resistivity 



where Ci and C2 depend 



is p(J) ~exp ^-^(log^) 

on various parameters such as temperature, magnetic in- 
duction, density and strength of point impurities as well 
as the strength of the defect plane. We conclude that a 
single relevant defect plane slows down the FL creep in 
comparison to the BG phase. 

There are interesting connections of some of the as- 
pects of our results to related two-dimensional classical or 
one-dimensional quantum models. A single planar defect 
in the Bragg glass phase resembles the presence of a sin- 
gle columnar defect in a FLL confined to a plan o^^'^^i^^ 
or a frozen impurity in a Luttinger liquid^i^. In all 
three cases the bulk phases on both sides of the defect 
are characterized by logarithmically diverging displace- 
ment correlations. The parameter g plays the role of a 
temperature in the 2D classical case and of the Luttinger 
liquid parameter in the ID quantum case. The periodic 
order seen around the defect plane has its counterpart in 
Friedel oscillations around an impurity in Luttinger liq- 
uids. Whereas in the ID (2D) case the relevance of an 
impurity is controlled by tuning the interaction strength 
(temperature), in the present case a change of g can be 
accomplished by changing the orientation of the defect 
with respect to the FLL. Transport properties of our sys- 
tem are however different from the ones in the related 
systems. 

In this article we also study the effect of a finite den- 
sity of randomly distributed parallel planar defects on the 
FLL at low temperatures with the magnetic field aligned 
parallel to the defects. We consider the case when the 
mean defect spacing is greater than FL spacing. Our re- 
sults may be directly applicable to a wide class of other 
systems with planar defects as a stack of membranes un- 
der tension, charge density waves^, domain walls in mag- 
nets and incommensurate systems^ since we consider a 
simplified model with an uniaxial displacement field per- 
pendicular to the defects. Wc find a new phase, which 
we call planar glass, that is characterized by (i) diverging 
shear (tilt) modulus that determines the energy cost for 
a shearing (tilting) of the FLL in the direction perpendic- 
ular to the planes; (ii) finite compressibility; (iii) sample 
to sample fiuctuations of the longitudinal magnetic sus- 
ceptibility; (iv) an exponential decay of positional corre- 
lations in the direction perpendicular to the defects and 
(v) a creep exponent ^ = 3/2 for creep in the direction 
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perpendicular to the defects for small currents J 0. 
The planar glass is different from the Bragg glass or the 
Bose glass phase and from the phase found for equally 
spaced defects^i. The planar glass is stable over a finite 
range of tipping angles of the applied magnetic field away 
from the direction parallel to the planar defects, i.e., it is 
characterized by a transverse Meissner effect. Similarly, 
the planar glass is characterized by a resistance against 
shear deformations that are perpendicular to the defects. 
Naturally, realistic samples contain both point and corre- 
lated disorder (as columnar and/or planar defects). We 
find that the planar glass is stable against both weak 
point and weak columnar disorder. The schematic phase 
diagram is shown in Fig. [TJ When we consider a vector 
displacement field, we additionally find that strong disor- 
der enforces arrays of dislocations in order to relax shear 
strain. They are located at the defects with a Burgers 
vector parallel to the defects. We argue that the main 
properties of the planar glass remain unchanged by dis- 
locations. 

The paper is organized as follows. In Sec. |lT] we in- 
troduce a model for interacting FLs that couple to weak 
point impurities and briefly review known results of this 
model. In Sec. |TTT] we consider a single defect plane as 
a perturbation to the Bragg glass phase and study the 
FL order using the a renormalization group analysis. A 
finite density of randomly distributed defects is explored 
in Sec. IIVI and the novel planar glass is characterized. 
Functional renormalization group equations are derived 
in c? = 6 — e dimensions. The response to tilt and shear 
deformations is discussed as well as sample to sample 
fluctuations of the longitudinal magnetic susceptibility. 
The positional correlation functions are computed and 
the stability of planar glass with respect to point and 
columnar disorder is studied. In Sec. |V] we consider the 
limit of strong planar defect potentials. In Sec. I VII we 
consider the FL dynamics for small currents by investi- 
gating FL creep in the presence of a single defect plane, 
both with and without point impurities, and in the pres- 
ence of a finite density of planes. Finally, in Sec. IVIl] we 
discuss a model with a vector displacement field. Tech- 
nical details and list of recurrent symbols are relegated 
to the Appendices. 



II. THE BRAGG GLASS PHASE 

In this section we summarize some known results on 
pinning effects due to randomly distributed point impu- 
rities for interacting FLs. We use elasticity theory to 
describe a dislocation free array of FLs (for a review see, 
e.g.. Blatter et al.— ). Undistorted FLs are exactly par- 
allel to the z-axis which we assume to be the direction 
of the applied magnetic field. The FLs form a triangular 
Abrikosov FLL in the xy-plane with a lattice constant a. 
In order to describe distortions of the FLs from the per- 
fect lattice positions R,y, we use a two-component vector 
displacement field u^{z). Since we are interested in the 



behavior on large length scales, it is appropriate to de- 
scribe the interacting FLs in terms of a continuum elas- 
tic approximation with a continuous displacement field 
u^(z) -> u(r). 
The Hamiltonian 

H^Ho + Hp (1) 

consists of the elastic energy of the FLL, Tio, and pinning 
energy of point impurities. Hp. The elastic energy of the 
dislocations free FLL reads 

(2) 

where qj^ = q^i + qy-y. = ft'Jj/?! and = 

S'ij — <li<lj/Qj_ are projectors onto the longitudinal and 
transversal modes, respectively, with propagators 

Gl^ = ciiqi + C44gf , (3) 

Gt^ = cgeqi + CMql- (4) 

In general, the compression (cn) and the tilt (C44) moduli 
are non-local on length scales smaller than the penetra- 
tion depth A but the shear modulus cqq is always local. 
However, the dispersion of cn and C44 on small length 
scales is negligible for the present problem, since we are 
interested in asymptotic properties at large length scales 
and small currents. Hence, in the following we introduce 
a cutoff in momentum space given by A « 2tt/X and ne- 
glect the non-locality of the elastic moduli. The elastic 
Hamiltonian of Eq. ^ can be obtained from symmetry 
arguments^. The ideal FLL is isotropic in xy— plane and 
has Dgfi symmetry group. 

The pinning energy of randomly distributed point im- 
purities is modeled by the coupling 

np = I rf3rp(r,u)yp(r), (5) 

of the local FL density p(r, u) = S{x — — u^{z)) 
to the pinning potential Vp(r), where x = ix,y). From 
this definition and the Poisson summation formula^i^ the 
density can be also written as 

p (r, u(r)) = po + Po I - V.u(r) + ^ e^G[x-u(r)] | 

(6) 

where po = 2/(\/3a^), and G is a vector of the recipro- 
cal lattice. Vp = —Vp ^^i^ ~ x,;)(5(z — Zi) represents 
the pining potential due to randomly distributed point 
impurities. The ^-functions are considered to have a fi- 
nite width of the order of the superconductor coherence 
length ^. For simplicity we subtract the average of the 
random potential and look at fluctuations around the av- 
erage value. The pinning potential then satisfies 

"M0=0, Vpir)Vpir') - n^mpvp^i^ - x')<5(z - z'). 

(7) 
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The strength of the disorder is characterized by v^riimp. 
Higher order correlations of the (unrenormahzed) pin- 
ning potential are nonzero but for weak disorder can be 
neglected. The restriction to two-point correlations of Vp 
leads to the same replica Hamiltonian one obtains when 
Vp would be Gaussian distributed. 

The model given by Eqs. ((I])-(I7]) has been stud- 
ied in detail using perturbation thcorjs^, Flory-typc 



argument; 
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a Gaussian variational ansata^iii^i and 



functional renormalization group^i^i^. The correlations 
of the FLL fluctuations change with length scale and are 
characterized by three different regimes: the Larkin or 
random force regime (RF), the random manifold regime 
(RM) and the Bragg glass (BG) phase. These regimes 
are distinguished by the scaling behavior of 



((u(r)-u(0))2 



cx r 



2C 



(8) 
de- 



which defines the roughness exponent C,. Here i 
notes a thermal and a disorder average. 

(i) In the Larkin regime^ the displacements are suf- 
ficiently small so that the FLs stay within one mini- 
mum of the disorder potential Vp(r) and perturbation 
theory can be applied. The effect of the disorder poten- 
tial on the FLL is properly described by a random force 
Fp(Ry,2:) = — VxVp(Ri,, 2). The roughness exponent 
is C-RF = (4 — 'i)/2, where d denotes dimension of the 
system, so that the positional correlation function 



Sg{v) = (e«Gu(r)g-»Gu(0)^ 



(9) 



decays exponentially fast in d = 3. The Fourier trans- 
form of Sg (r) is the structure factor which can be directly 
measured in diffraction experiments. The Larkin lengths 
L| and L| are defined as the crossover length scales 

where the conditions ((u(0,z = L|) — u(0))2) cx and 

((u(|x| ^ L|,0) - u(0))2) cx ^2 are satisfied. This leads 
to 
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(j)o£,^ 044066 

Bv'^Uimp 1 + K 

A e6 ^l/2„ 3/2 
<Pot, C44 '^66 

Bvjn^rnp 1 + K^/^ 



(10) 



where (po = hc/{2e) ~ 2.07 10^ G cm is the flux quan- 
tum and K = cee/cii. The length increases with de- 
creasing disorder strength. An increase in the magnetic 
induction B effectively increases the disorder strength so 
that L^ shrinks. 

(ii) On scales greater than the Larkin length, a de- 
scription in terms of random forces become inapplicable 
and the RM regime applies. In the RM regime FLs ex- 
plore many minima of the disorder potential but the typ- 
ical displacement is still smaller than the FL spacing a. 
Hence the FLs do not compete with neighboring FL for 
identical pinning centers. A Flory-type argument^i^S 
yields the roughness exponent C-RM = (4 — d)/6 but in 
Rcf. jSj] is was shown within an e = 4 — d expansion that 



Crm depends on the ratio k = Cqq/ch and varies be- 
tween 0.1737e and 0.1763e. The positional order decays 
according to a stretched exponential, 



Sair) ~ exp 



(11) 



(iii) On length scales larger than the positional corre- 
lation length La the RM regime becomes inapplicable. 
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'/^'^ {a/S,)^^^'^" is defined as the scale at which 
the mean square displacement of FLs is of the order a. 
Therefore, it is crucial to keep the periodicity u u-f Ri, 
of the interaction between FLs and point disorder—. This 
leads to a much slower logarithmic increase (Cbg = 0) 
of the clastic displacement of the FLs than in the Larkin 
and RM regime. It was showrt^'^'^ i^i^i^'^ that thermal fluc- 
tuations are irrelevant and that the pinned FLL exhibits 
a power law decay of positional correlations. 



5g(x,0) - |x| 



(12) 



where j/g = v{G/GoY and Go = 47r/(\/3a). This re- 
sult resembles the correlations in pure 2D crystals at fl- 
nite temperatures. A functional renormalization group 
analysis in 0? = 4 — e dimensions yields a non-universal 
exponent 77 that varies with the elastic constants of the 
FLL^i^. Extrapolating to d = 3, one finds only a very 
weak variation with 1.143 < 77 < 1.1592i£. Despite of the 
glassy nature of the phase algebraically divergent Bragg 
peaks still exist which motivated the name Bragg glass^ii. 
The existence of the Bragg glass phase has been experi- 
mentally confirmediSii^. 

After this summary of the scaling regimes, we briefly 
review the replica theory for the Hamiltonian of Eq. Ilj . 
Using the replica method, we average over point impuri- 
ties (see Appendix A) and obtain the replica Hamiltonian 

n n ^ 

n-p = Yl ^oK) - ^ E J i?pK(r) - u''(r)]. 

(13) 



Q,/3=l ■ 



Rp{u) = {vppafn,mp ^""^"k-^ (G), (14) 

where (5j-i(G) is the delta function smeared out over a 
region of size The correlation functions Ct and Co 
that describe thermal fluctuations and disorder induced 
fluctuations can be written as 



GT = (S(q)S(-q))-(ii(q))(S(-q)) 
= (27r)'^r{gi(q)+^T(q)} 



Cd = (5(q))(S(-q)) = (2^)'^A(q){ei(q) +e|(q)} , 

(15) 

where the last equation defines A(q), which we shell ob- 
tain in harmonic approximation below. 
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In the next section we study the interplay between 
point impurities and a planar defect. This is a difh- 
cult problem since we have to deal with two nonlinear 
terms. We consider the planar defect as a perturbation 
to the BG fixed point and examine the stability of the 
BG phase. Also we explore the effects of the defect on the 
order of the FLL. In the following we will use an effective 
quadratic Hamiltonian that reproduces the displacement 
correlations of Eq. ^ of the full nonlinear disordered 
model Eq. p^ . A systematic analysis must be based on 
e = A — d expansion, and a functional renormalization 
group analysis shows that displacements obey Gaussian 
statistics to lowest order in e^. It should be noted that 
the effective Hamiltonian does not capture all physics, in 
particular, it cannot describe correctly the FL dynam- 
ics since it cannot reproduce the energy barriers for FL 
motion''^. An effective quadratic Hamiltonian has been 
also used for a model with an uniaxial displacement to 
study a dislocation mediated transition of the FLL4^. 

The effective quadratic replica Hamiltonian in d di- 
mensions reads^i^ 



parallel to the z axis, we will calculate all the integrals 
in d = 3 if not stated otherwise. If the numerical values 
of rj and Crm are important for our conclusions, we will 
comment on a possible influence that the use of results 
found by Emig et ali^i^ can have. 



III. SINGLE DEFECT 

In this section the influence of a single planar defect 
on the Bragg glass order of the FLL is studied. In some 
parts of this section, when examining FLs density os- 
cillations around the defect, we will study the isotropic 
limit with cn = C44 = cee = c in order to focus on the 
important physics. In this limit the propagators read 
CJ^^(q) = tjy^(q) = cq^. By this assumption, only the 
weak dependence of 77 and C,b,m on the elastic constants 
is ignored. 



1 

- y 

Q,/3 = l 



(2^) 



-d 



(16) 



where 



QJ:M) =<5„,^fel(q)Pi +C;^^(q)P5 



.^1 

T 



A(q) 
T 



1. 



(17) 



It yields the correlation functions of Eq. p^ . where A(q) 
describes the behavior of A = — i9^^i?p(0) = —9^ i?p(0) 
on different length scales. Using a functional renormal- 
ization group in d = 4 — e dimensions, it has been shown 
that to lowest order in e^i^ 



A(q) 



1, 



(18) 



The function A(q) reaches the fixed point form 
g'^A*(K)c44C66a^ in the BG phase, where A*(k) ^ e/(l + 
k) depends only on elastic constants but not on the dis- 
order strength. We note that Emig et al^i^ have ob- 
tained their results by calculating the integrals, needed 
for the RG equations, systematically for c? = 4 with a 
two-dimensional vector z. The results are then extended 
to three dimensions by setting e = 1 in A*(k). This 
approach does not influence the main physics (like the 
logarithmic roughness of FL in the BG phase), but may 
influence the dependence of exponents 77 and Ci?J\/ on the 
elastic constants. In this way the dimensionality of z 
" axis" and the contribution of the term cac^^ in the prop- 
agators are more weighted than the other axes and other 
terms ~ C66,cii, respectively. In the following, in order 
to not overestimate the effect of a planar defect that is 



A. Model 

The pinning energy of a planar defect can be written 
in the form 

'Hd = j d'rp(r, n)VD (r • - (5) , (19) 

where VD(r • — 5) is the potential of the defect plane, 
no and 5 denote the unit vector perpendicular to the 
defect plane and its distance (along n^i) from the origin 
of the coordinate system, respectively. The Bragg glass 
order that we are interested in is dominated by disor- 
der fluctuations on large length scales where microscopic 
details become irrelevant. Therefore we may approxi- 
mate the defect potential by a smeared out (5-function, 
Vd{x) w —v5^{x). Since the superconducting order is re- 
duced in the defect plane, it is plausible to assume w > 
(for more details, see Section IX of Blatter et al.— ). When 
we assume that FLs gain condensation energy when they 
overlap with the defect plane, a rough estimate for the 
defect strength is v k, H^^^ with He the thermodynamic 
critical field. 

In order to integrate over the delta function of the de- 
fect potential, it is convenient to introduce an explicit 
parametrization for the position vector r^i of the defect 
plane which obeys ru-no = S. With the parametrization 

ro = (x_d,z_d) +(5n_D, ZD^tcos/3 (20) 
X£) = (s sina — t COS a sin/3, s cosa -|- t sina sin/3) 
nu = (cos /3cos a, — cos/3sinQ!, sin /?) 

we introduce in-plane coordinates s, t, and the two angles 
a and P which determine the rotation of the plane with 
respect to the y- and z-axis, respectively (see Fig. 
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The defect energy now reads 

TYd = vpo J dt ds dr±S^{r± — ^)|Vxu(t, s, r±) 

_ y ^ ^iG[r±no+xo-u(t,s,r±)])^ {21) 

where r± = r • njj. Since the displacement field u varies 
slowly on the scale of the FLL constant a, the integrals 
over s and t vanish for all G with the exception of those 
for which the oscillatory factor e^'^^'^ is unity (for all s, 
t). This condition can be satisfied only if sin (3 — 0, i.e., if 
the defect plane is parallel to the applied magnetic field. 
There remains a second condition for the angle a which 
results from the constraint that G — mbi + nb2, with 
integer m and n, has to be perpendicular to x^i. Express- 
ing the defect plane (for sin /3 = 0) as xd — (ciai — C2a2)s 
where a^bj = inSij, one sees that the second condition 
is equivalent to the condition m/n = c^jcx. Hence if 
ci/c2 is irrational, the effect of the defect plane is al- 
ways averaged to zero. On the other hand, for ratio- 
nal cijcx we may choose mTj^nj;) to be the smallest co- 
prime pair with C2/C1 = m]j/n£,. Then too, nj) are 
the Miller indices of the defect plane and only those G 
which are integer multiples of G^ = rn^bi -f noi>2 con- 
tribute in Eq. ((2T|) . In the following, we will concentrate 
on the contribution from these G-vectors only. The FLL 
planes (of the ideal lattice) that are parallel to a defect 
plane with Miller indices to^, ud have a separation of 
—a/ \/m'jj + monD + «|) and hence Gd — 27r/£. 



2 ^1 y + moriD + njj and hence Gd 
For the defect plane aligned to the magnetic field we 
take the a;-axis to be perpendicular to the defect (i.e. a = 
(3 = 0) and hence the defect Hamiltonian becomes 



/L^/ S J 
dydz ^^V ^u(r d) " ^ 2 cos [fcG'D(<5 - Uj;)] |, 
fe>o 

(22) 



where = {S,y,z). Here Ux denotes the component of 
the displacement field that is perpendicular to the defect 
plane and [x]q is the integer number that is closest to x. 



B. Renormalization group analysis 

In this subsection we discuss the influence of the planar 
defect on the stability of the BG phase using a renormal- 
ization group (RG) analysis. We employ a sharp-cutoff 
scheme by integrating out the displacement field (q) , 
with wave vectors q in an infinitesimal momentum shell 
below the cutoff A > |q| > A/b = Ae~' and subsequently 
rescale lengths and momenta according to 



qfe 



(23) 
(24) 



We split the displacement field into weakly varying modes 
M^(r) and strongly varying modes u^(r) that include 




FIG. 2: Vectors of the triangular flux line lattice (ai, a2) and 
of its reciprocal lattice (bi, b2), and the angles a, /3 that 
define the orientation of the defect plane. 



Fourier components out of and in the momentum shell, 
respectively. We choose to not rescale the field u'(r') = 
u^(r) which implies a rescaling of its Fourier transform, 
ii'(q') = S<(q)/63. 

The defect plane is considered as a perturbation to the 
Hamiltonian of Eq. The gradient term of Eq. ((^ 

scales ~ L if the defect size ^ L^. Since the elastic 
energy Eq. ([2]) scales in the same way, the gradient term 
is a marginal perturbation. It can be also eliminated by 
the transformation u'^{r) = ''^x{^) + sgn.{x — S), where 
sgn(O) = 0. This transformation does neither change the 
terms ~ cge, C44 of Eq. ([2|) nor the pinning energy due 
to point impurities in Eq. (fT3|) since all replica fields are 
transformed in the same way. The gradient term of the 
defect pinning energy tends to increase the FL density at 
the defect as can be seen from the transformation above. 

In order to account for different renormalization of the 
harmonic components of the defect pinning energy, we 
introduce the variables Vk for the strengths of the har- 
monics of order k. A cumulant expansion yields to first 
order in v the renormalization 



Vkjl) _ v_ 
T{1) ^ T 



e2'(cos [kGoK'^iro)]) 



T 



\{kGD?{[uZ-'^(rD)?) 



(25) 



where the factor e is due to a rescaling of lengths. 
([u"^>(r£i)]^) is obtained at the BG fixed point to linear 
order in I. Due to the irrelevance of thermal fluctuations, 
we have neglected contributions that come from the ther- 
mal part of the propagator of Eq. (|17|) . We have chosen to 
rescale temperature instead of elastic constants in order 
to organize the RG analysis of the zero temperature BG 
fixed point. It is important to note that Eq. ((25| holds 
only on length scales larger than the positional correla- 
tion length La- 
in the random manifold regime (e*'^"('")) decays with 
the system size as a stretched exponential and the effect 
of the defect plane is reduced by disorder fluctuations on 
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FIG. 3: Two possible orientations of defect planes relative to 
the flux line lattice, corresponding to g — 77/2 (dotted lines) 
and 3r;/2 (dashed lines). 



intermediate length seales. Hcnee, the renormalized and 
rescaled value of the defect strength is reduced to 



(26) 



on the scale L = La, where C is a positive constant. This 
value is the initial defect strength Vk{l = 0) to be used in 
Eq. (HID. 

The RG flow equations in the Bragg glass regime now 
read 



dl 

dVk ,2 ^ 



(27) 
(28) 



Hence Vi is a relevant perturbation provided g < 1, i.e., 
if 



0.66 a. 



T]{m.^ + moriD + rijj) < 2 or ^ > y -3- 




(29) 

which is compatible only with £ = \/3a/2 « 0.87a. A 
relevant defect plane must be oriented parallel to one 
of the three main crystallographic planes of the FLL 
(i.e. cos 2/3 = cos 6a = 1). When £ increases {g de- 
creases) more FLs can gain energy from the defect plane 
and hence render it more relevant. 

Emig et al^i^ have calculated rj in a one-loop functional 
RG expansion in 4 — e dimensions. Higher loops as well 
as the fact that all integrals are evaluated in d = 4 with z 
being a two dimensional vector, may influence the actual 
numerical value of the coefficient rj in d = 3. However, 
we argue that this higher order correction does not affect 
our conclusion that a single defect is relevant only if it 
is parallel to the main crystallographic planes, since g = 



(77/2, 377/2, 777/2, 
Fig.©. 



7-1, )/2 can change only in finite steps 
. ) when rotating the defect plane (see 



C. Effective Hamiltonian 

In this section we discuss whether higher order cumu- 
lants in v can lead to a renormalization of the parameter 
g and hence can influence the condition for the relevance 
of a defect that was derived in the previous section. The 
renormalization described by Eq. (|28p does not occur in 
the bulk but on the defect plane. Hence it is possible 
to develop an effective theory that is defined on the de- 
fect plane only. Since the defect couples only to the dis- 
placement Ux on the defect plane, we integrate out 
outside the defect and Uy across the entire sample. This 
integration is facilitated by employing the effective Gaus- 
sian theory for the BG phase of the previous Section. At 
T = we are interested in the ground state and hence we 
solve the Euler-Lagrange equation for u(r) with the con- 
dition Ma;(r£)) = ip{rD) at the defect plane, where firo) 
is an arbitrary function. An equivalent functional inte- 
gral approach is presented in Appendix [B] The effective 
replica Hamiltonian on the defect plane reads 

^e// = - EE ^^'^-Po / dvD cos {kGoiS - p"{rD)]} 

a k>0 •' 

(30) 



a, 13 



where q is the in-plane momentum. 



(S^(x,q)u^(x,q')) =r(27r)'*<5(q + q')Q "^(q) (31) 

and Vk are renormalized parameters on the scale of the 
positional correlation length. In order to avoid technical 
complications, we consider the limit of isotropic elasticity. 
In line with an e expansion, we evaluate the integrals in 
d = 4 and then set e = 1 in the expression for fixed 
point value A*(k). This approach does not affect our 
conclusion and leads to a clearer result. On scales larger 
than La we get 



T 
T 



cq 



cnA 



BG 



T 



[qcf 



(32) 



where c is the elastic constant and A^g ~ c^a^ AA*(1). 
The same procedure can be performed also in the RF and 
RM regimes using the corresponding quadratic Hamil- 
tonian in d dimensions. The effective Hamiltonian in 
d — 1 dimensions has a long ranged elasticity (term ^ q 
in the limit Abg 0) that results from the local bulk 
elasticity. A RG analysis of the effective Hamiltonian 
of Eq. pO|) shows that neither Abg nor the elastic con- 
stants are renormalized, and hence g is not renormalized. 
From this we conclude that a weak defect is a relevant 
perturbation only for g < 1. 
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D. Density oscillations 



In the limit T ^ this result can be expressed as 



Next, we study the order of the FLs next to a de- 
fect plane. We consider separately the case of a rele- 
vant and an irrelevant planar defect plane. For simplic- 
ity, we assume isotropic elasticity and choose to place 
the origin of the coordinate system on the defect plane, 
i.e., wet set S = 0. In the absence of a planar de- 
fect, PL den sity fluctuations due to point impurities obey 
{p{r) — po) = 0. The defect plane pins FLs and yields a 
long-ranged restoration of the translational order param- 
eter e*'^"('"). We find Friedel-like oscillations of the FL 
density with an amplitude that decays as a power law 
with an exponent that depends on if the defect is rele- 
vant or irrelevant in the RG sense. 



1. Irrelevant defect 

First, we consider an irrelevant defect parallel to the 
magnetic field. The irrelevance of the defect potential 
for g > 1 allows us to compute the thermal and disor- 
der average of FL density perturbatively in the defect 
strength. 



(5p(r,u(r))) = lim TT / Pu" 5p(r,u''(r)) e' 

n— >0 J 



(33) 



where Sp{r,u) = p(r, u) — po, 7 is an arbitrary replica 
index and 

=HJ^ + ^7^i,(u"). (34) 

a 

To lowest order in the defect strength we get 



{6p{r, u(r))) = lim {6p{r, u''(r))) ^ . (35) 

Even if the defect is irrelevant in the RG sense, it breaks 
the translational symmetry perpendicular to the defect 
and hence modifies the FL density locally. To cor- 
rectly describe this effect we need to compute the average 
change in the FL density to first order in v. We find (see 
App.p 

n 

(<5p(r,u(r))> =-/3 lim V(<5p(r, u''(r))Hz,(u")). (36) 

71—^0 ^ ^ 



(Mr,u(r))) « ^-ld^,osiGo\x\) 

c{2g-l) \\x\ 



2g~l 



(37) 



which becomes exact for large x. The result captures the 
large length scale behavior for L > La- Here vi denotes 
the effective defect strength on the scale La, cf. Eq. (|26p . 
and is the normal distance from the defect plane. 
There are additional contributions to Eq. (P7|) coming 
from the higher harmonics in TY/j. They are less impor- 
tant since they are proportional to w^, on the scale La 
and they decay as |a;|~^'^ ^"""^ with k > 2. Although the 
defect is irrelevant in RG sense, it leads to Friedel-like 
oscillations in the density. 

If the defect plane is not parallel to the applied mag- 
netic field, Friedel oscillations occur as well. However, 
the amplitude is exponentially suppressed. The ampli- 
tude of density oscillations with reciprocal lattice vectors 
G = G(cos a, — sin a) (for a definition of the angles a and 
f3 see Fig. [2]) decays beyond the distance I/(G|sin/3|) 
from the defect plane. Similar physics occur in classical 
2-dimensional systems with a columnar defec t^^'^^ . 



2. Relevant defect 

The strength of a relevant defect grows under renor- 
malization relative to the elastic and the impurity energy. 
On the scale 



vLn 



(38) 



the energies become of the same order and perturbation 
theory breaks down. On larger scales, the defect po- 
tential can be described effectively through the bound- 
ary condition Ux{x = 0,y,z) = for the displacement 
field at the defect plane. With this constraint the sys- 
tem gains maximal energy from the defect and the com- 
plete energy of the system is minimized. First, we cal- 
culate displacement correlations Qpin,ijix,x';r\\ ~ ''y) ~ 

{ui{r)uj{r')) with the above boundary condition at 
the defect and r ~ (a;,r||). We find in momentum space 
(see Appendix [D|) 



iix,x;q) lirn ^"(0, q) - ^ g^"(|a;|, q)Q„_^^(q)g^](|x|, -q) 



(39) 



9 



~a,P 1 

where q is the m-plane momentum and Q is the inverse of Qajj given by Eq. pTj) . i,j stand for x,y and Q is 
given by Eq. (PT|) . It can be shown that the displacement correlations on scales larger than Ly are given by 



{ui{r)uj{v)) 



lim T 



g"(0,O)-(i.i)(i.j)(?ii(2|x|,O) 



(40) 



Using this result, we obtain for the average change in the 
FL density 

(Mr,u(r)))=2po5]cos(mGDx)^j . (41) 

m>0 ^ ' ' ^ 

These oscillations resemble Friedel oscillations which 
can be also found in Luttinger liquids with an isolated 
impurity^ or in classical 2-dimensional systems with a 
columnar defec t^^i^^ . 

The amplitude of the Friedel oscillations decays as a 
power law with an exponent g and 2g — 1 for a relevant 
and an irrelevant defect, respectively. For an irrelevant 
defect the amplitude decays more rapidly than for a rele- 
vant defect. In the absence of point impurities the defect 
is always relevant in the RG sense, and the amplitude of 
the Friedel oscillations remains finite for |.t| — > oo. 

IV. FINITE DENSITY OF WEAK DEFECTS 
A. Model 



defects and point impurities. The defects tend to local- 
ize the FLs and hence favor order along the defect planes 
while point impurities promote FL wandering. 

The Hamiltonian reads 

n = no + j d\ [Vp (r) + 1/d (r)] p(r , u) , (42) 

where TLq is the elastic Hamiltonian of Eq. ([2]) and Vp 
is the pinning potential resulting from point impurities, 
see Eq. The defect pinning potential is Vd{y) = 

—v{^^ S^{x — Xi) — 1/£d} where we assumed that the 
defect planes are parallel to the j/z-plane and £d is a 
mean defect spacing. The i5-function is assumed to have a 
finite width of the order of the superconductor coherence 
length ^. The defect potential is uncorrelated along the 
X-axis, 



VD(ri)Fz?(r2) = ^S^{xi - X2). (43) 



In this section wc consider a finite density of parallel Wc discuss the case where the gap between two defect 



planar defects with random position. We assume that 
defects extend along the entire sample and arc aligned 
parallel to the applied magnetic field. There is a com- 
petition between the two random potentials from planar 



planes typically contains many FLs, i.e., ^ <C ^ ^ ^_d. 
Note that orientation of the defects is otherwise arbitrary. 

After averaging over the defect positions, the replica 
Hamiltonian for the defects reads 



J 



D 



2T£n 



/"d3rid3r2^<54(x-i-X2)| -2Vxu"(ri) 



G 



Vxu"(ri)VxU^(r2) + ^ 



giGi [xi-u° (n)] giG2 [x2-u'' (r2)] 



(44) 



Gi,G2 



where x = {x,y). The defects are assumed to be sufficiently weak so that terms of the order v'^ /io and higher can be 
neglected. The first term does not contribute to due to the oscillatory factor e^^^^ and the third term contributes 
only for reciprocal vectors perpendicular to the defects satisfying Gi = — G2 = uGd with integer n. Introducing 
the relative coordinate Xr = xi ~ X2 and taking into account that S^{xr) is finite for \xr\ < v/e approximate the 
displacement field as Ux{x2 + Xr,yi, zi) w Ux{x2, yi, zi). This approximation is justified since the displacement field 
varies slowly over the FL spacing. Then Eq. (|44p can be written as 



^ " 2T 



{(TVxu"(a;i,?/i,zi)Vxu''(xi, 2/2,^2) + Rd [u"{xi,yi, zi) 



uf(xi,y2,Z2)]} , (45) 



where and wc defined = J dx and ct = {vp^Y / Ijj. After aver- 

Rd{u.) = aJ2 h-' (nGz?)e"^^-"^ (46) 
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aging over point impurities, the complete replica Hamil- 
tonian is TT ^Wp+TTiy, where is given by Eq. (fTS]) . 

The first term in Eq. ([^5]) comes from the coupling 
of the defect potential to the slowly varying part of the 
FL density ^ VxU. This term does not contribute to the 
glassy properties of the system, since it can be eliminated 
by a simple transformation^)^. (For a more detailed 
discussion of this term see below.) The remaining part 
of the replica pinning energy TL" is invariant under the 
transformations 



<(r)-.<(r) + /.(x), 
<(r)^^.^(r) + /,(r), 



(47) 
(48) 



where fx{x) and /y(r) are arbitrary functions. Eq. ((47|) 
represents an approximate symmetry if the defect poten- 
tial has a finite width. However, with increasing length 
scale, deviations from the symmetry become less impor- 
tant. These symmetries show that the elastic coefficient 
cii is not rcnormalized. Not renormalized arc also the 
elastic moduli which determine the energy cost for tilting 
the FLs only in the y direction (i.e. the term c^AidzUyY) 
and for changing only the displacement Uy along the x- 
axis (i.e. the term CQ^idxUyY). These symmetries are 
commonly denoted as statistical tilt symmetr}*^^. How- 
ever, the defects are an important source of anisotropy 
and other elastic properties of the FLL will be affected. 
For example, the energy cost for FL tilting as well as FLL 
shearing parallel and perpendicular to the planes will dif- 
fer considerably. Also, due to the defect planes the sys- 
tem is not invariant under arbitrary rigid rotations of the 
FLL around z axis and rotational modes will appear in 
the elastic Hamiltonian under renormalizatior4^. Note 
that for the FLL with point disorder only, none of the 
elastic constants will be renormalized since the disorder 
correlation function Rp{\i) is invariant under the more 
general transformation u" (r) u" (r) -|- f (r) . 

Planar defects in the form of twin boundaries that are 
perpendicular to the copper oxide planes very often ap- 
pear in YBazCusOy-o: (YBCO)i^ii^ii^. YBCO is a high 
temperature superconductor and within high accuracy 
it is uniaxially anisotropio^. YBCO can be reasonably 
well described within a continuum anisotropic model, 
while for more strongly layered superconductors a dif- 
ferent description is needed. The elastic description for 
anisotropic superconductors can be found in the review 
article by Blatter et al.— . The number of independent 
elastic moduli increases with respect to the isotropic case 
that we discussed in Sec. [Hi However, if the magnetic 
field is applied perpendicular to the copper oxide planes, 
the model given by Eq. ^ as well as the considerations 
in the following sections are directly applicable also for 
YBCO. 



consider both the planar defects and the point impurities 
as a perturbation to the ideal Shubnikov phase. Notice 
that Eq. (j^ depends only the displacement field u^- 
Since we focus below on the effect of the defect planes, 
it seems to be justified to start from a simplified model 
in which only the displacement Ux = u oi the FLs per- 
pendicular to the defect planes is considered. This model 
describes also a wide class of other systems which exhibit 
regular lattices of domain walls like magnets, charge den- 
sity waves^ and incommensurate systems^^. 

In the absence of the defect planes point impurities 
are relevant below four dimensions. We employ here an 
Imry-Ma-type argument^ in combination with pertur- 
bation theory to see the effect of randomly distributed 
point impurities on the FLL. When the initially ordered 
FLL is distorted in a volume L'^ by u ~ ^, the typical en- 
ergy gain is of the order ~ {—Rp{0)L'^y^^ compared to 
the clastic energy loss ~ L'^~^. For d < dp ^ 4 and suf- 
ficiently large L:^ L^r^ {-R'^{0)f'^'^^^^ > (. the point 
disorder wins and the FLL becomes distorted. A more 
detailed study shows that in this case the FLL exhibits a 
phase with quasi long range order which is the previously 
discussed Bragg glass phase (see Sec. In this phase 
the positional correlation function Sc, [see Ec^. shows 
a power law decay. 

Next we consider the Imry-Ma argument for planar 
defects in a volume L'^~^LzLy without point impurities. 
The energy gain is of the order {-R'ij{0)eLi-'^f''^LzLy. 
The elastic energy loss is ciiLzLyL'^~'^£'^ since distor- 
tions are aligned parallel to the defects. For Lx ^ Lp) ~ 
(— cf]^^^/i?^(0))^/'^~'^' the pinning energy gain wins and 
the FLL becomes disordered in the direction perpendic- 
ular to the defects. Lp) is the so-called Larkin length 
for the defects. The critical dimension above which weak 
planar defects are irrelevant is dp — 6. 

For an RG approach is convenient to consider a gen- 
eralization of our model to d dimensions. The defects 
remain two-dimensional with d — 2 transverse directions, 
while the displacement field remains uniaxial. In the fol- 



lowing, we use a functional RG approach 



dimensions. We follow closely a related approach for 
columnar disorde r^^i^° but do not rcscale the renormal- 
ized quantities so that they correspond to the effective 
parameters measured on the scale Lx- Thermal fluctu- 
ations and point disorder are irrelevant for e < 4 and 
e < 2, respectively. Hence we can assume directly T = 
and Rp = 0. To lowest order in e the RG flow equation 
read 



dine,, ^ KdR'j!,' {0)L% 
din La; 
dRD{u) 



-11 



d in L-r 



KdR!i,{u)L% 



4,6 



[R'Uu)~2R'i,{ 



(49) 
(50) 



B. Functional renormalization group approach 

In the previous section we treated a single defect plane 
as a perturbation to the Bragg glass fixed point. Now, we 



where Kd = Sd-2/ {2i^) ^i^d Sd denotes the surface of 
the d-dimensional unit sphere. 

In a static situation the displacement field is indepen- 
dent on y and z since defects distort FLL planes that 
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are parallel to the yz-plane on the whole. Since, by as- 
sumption, other sources of fluctuations are not present we 
can perform the integration over y and z in the Hamilto- 
nian of Eq. (|^5|l and obtain an effective d—2 dimensional 
Hamiltonian that describes the interaction of FLL planes 
with defect planes with d—2 dimensional random posi- 
tions. This explains why the flow equation for has 
the form as the one for the point disorder correlator Rp 
in d — 2 dimensions^. However, an important difference 
between the d dimensional FLL with defects and the FLL 
planes with point impurities in d — 2 dimensions is the 
renormalization of elastic constants C44 and cge in the 
former model. 

For Lx L]j, Rjj (0) increases and at ~ Ljj Rjy{u) 
develops a cusp at u = 0. The cusp signals the appear- 
ance of metastable states the energy of which is very close 
to the ground state energy but which may be far apart in 
conflguration space. The cusp results in diverging elas- 
tic constants C44 and cge and in a change of the sign of 
i?""(0+) from positive to negative. If there is a small but 
finite tilt or shear of FLs, i?^ (0) has to be replaced by 
R""{0+) in Eq. and on length scales > Ld the 

elastic constants decrease since i?""(0+) is then negative. 
Importantly, a new term of the form 
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This result is in qualitative agreement with the estimate 
we obtained from scaling arguments since Rjj (0, X)£'^ ~ 
i?;(0,A). 

The fixed point function has for < u < £ the form^ 



R*jj {u,L,j:) 



6K, 



""2j ^12 



(55) 



jd-2 



xdydz\Y.y{dyu)y + 'E;:{dzu)z\ (51) 



and it has to be periodically continued in u with period £. 
Note that if we would consider rescaled quantities then 
in Eq. ([55]) L~'^ would be replaced by A*^. In the case 
of a finite tilt (shear) for Lx > Ljj the elastic constant 
C44 (cge) decreases as Cii{Lx) ~ {Lx/ Ld)~'^^^ ■ A rough 
estimate for the saturation value for the interface tension 
is e£'^^JciiCii{\)/LD. 

As has been pointed out by Fedorenko^, the RG equa- 
tions for the elastic constants C44, cge and the interface 
tensions resemble those of the friction and driving force 
for the depinning transition of the FLL in the presence of 
point impurities^. The role of velocity is here played by 
tilt or shear and the elastic constants diverge at zero tilt 
and zero shear as the friction diverges in the static case. 
Also, as a threshold force exists for the depinning transi- 
tion, the interface tension ^^(.y) determines the threshold 
force for tilting (shearing) the FLL, as we will see in the 
next subsection. 



is generated in the Hamiltonian. '^z(y) has the meaning 
of a interface tension of a domain wall perpendicular to z 
(y) axis, across which the displacement field changes by 
£. They dominate the elastic energy for small u and are 
renormalized according to 



-1/2 



dS„ 



-1/2 



dS„ 



-66 



din Lx 



din Lx 



^Ro{0+)Ll: 



93/2 



(52) 



S2 and Yiy satisfy the relation Ti^/'Siy = ^/cil/^/cge. No- 
tice that C44 and cqq are renormalized in the same way 
such that their ratio remains constant under the RG flow. 

Next we estimate the Larkin length Lo from the flow 
equations. The function Rd{u) is even and as long as it 
is analytic, all odd derivatives at u = vanish. Assuming 
analyticity, the flow equation for i?^ (0) reads 



d In L r 



-KdL%[Rj, (0)]^ 



-^11 



The 

Lx 



solution on the 
is given by Rr 



(53) 



scale 



length 

D (0,^4^ 

(0, A) [1 - i?^' (0, \)^Kd {L% - /cl^e] 

A is the penetration depth and has a role of the small 
length scale cutoff. This shows that Rd{0) diverges at 



where 



L 



D 



AKdR'oiOA) 



(54) 



C. Properties of the planar glass 

In this subsection we summarize the properties of the 
new phase that is described by the fixed point of the 
functional RG of the previous section and in the following 
is called "planar glass". We examine the response of 
the system to FL tilting, to a change in the longitudinal 
magnetic field and discuss the order of the FLL. We show 
that the new phase and its properties are robust against 
weak point impurities in d = 3 and d = 4. 

When one changes the direction of the applied mag- 
netic field by Hx^i-, the Hamiltonian changes by 



An 



d^r HxdzU. 



(56) 



To tilt the FLs with respect to the z-axis, Hx has to 
overcome the interface energy ~ which results in a 
threshold field 



Hx c — 27r\/3' 



E^ 



(57) 



below which the FLs remain locked parallel to the planes. 
This is the transverse Meissner effect: a weak transverse 
magnetic field Hx is screened from the sample and C44 
is infinite. Only for Hx > Hx^c the average tilt of the 
FLs becomes non-zero and C44 is finite. In this way, by 
measuring the threshold field, can be measured. 

Moreover, there is a resistance against shear of the 
FLL. The shear deformation dyUx is non zero (and cgg is 
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finite) only if the shear stress (Jxy is larger than a critical 



value axy 



Otherwise cge is infinite. The diver- 



gence of C66 is a new property that does not appear in the 
Bose glass which, however, does also show a transverse 
Meissner effect. 

An infinitesimal change in the longitudinal magnetic 
field SHzZ changes the Hamiltonian by 



in ' 



(58) 



and allows to measure the longitudinal magnetic suscep- 
tibility X = <PoPad{dxu) / d5Hz. The disorder averaged 
susceptibility is 



X 



{<l>aPof 



(59) 



as shown in Appendix [E] It is independent of disorder 
as a result of the statistical tilt symmetry^. The glassy 
properties of the system can most easily be seen by the 
sample to sample fluctuations of the magnetic suscepti- 
bility. Perturbation theory yields (see Appendix IE)1 



X - X 
x' 



Lx_ 
Ld 



(60) 



i.e. the sample to sample fluctuations of the susceptibility 
grow with the scale ^ Lu, d < 6. We cannot expect 
that this result is quantitatively correct for large , but 
qualitatively it demonstrates the relevance of defects and 
it provides a signature of a glassy phased. Although 
we were not able to prove it, x^/x^ ~ 1 will most likely 
approach a finite universal value for 3> Ld in d < 6. 

The positional correlation function is obtained to first 
order in perturbation theory, combined with a functional 
RG analysis for 6 > d > 4. It reads 



-5*00 2) ^ \A 



(61) 



where ryu = (7r/3)^(6 — d). A detailed derivation of this 
result is presented in Appendix |F1 In order to study 
the behavior in d < 4, we have to reconsider the first 
term of Eq. (pSj) . It results from the coupling of the de- 
fect potential to the slowly varying part of the FL den- 
sity (^ dxU). By taking into account that at a scale 
the displacement field behaves as w ~ L£, we find that 
the (T-term scales as ^ When we com- 

pare the latter term to the squared elastic Hamiltonian 
^ LyL^Lx^''' that describes the cost of deviations 

of u in the x directions only (since the FLs are completely 
ordered parallel to the defects on sufficiently large scales 
in the absence of point impurities), we find that the a- 
term becomes relevant if d — 4 -I- 2C < 0. For logarithmic 
roughness (C = 0) it is relevant for d < 4. Since the 
other part of the defect pinning energy ^ Rd{u) scales 
in the same way as the elastic energy, the cr-term is the 
dominant part of the pinning energy and determines the 
FL roughness. 



First, we consider the case without point impurities 
and then treat them perturbatively. Applying a Flory- 
type argument^, i.e., assuming that the elastic energy 
and the dominant part of the defect energy scale in the 
same way, wc find, following the discussion above, that in 
d = 3 the roughness exponent is C = 1/2. More detailed 
calculations^!^ confirm our result, leading to 



,(x,?/,z) - e 



-kl/?c 



(62) 



where ~ i_D- Note that there is a shift of dimen- 
sion d ^ d -|- 2 between the model studied in Refi^i^ 
to our model since the FLs are ordered in the yz-plane. 
In Ref.i^>^ a related one-dimensional system with point 
impurities at zero temperature is studied. There is a non- 
trivial renormalization of a coming from the defect poten- 
tial that couples to the periodic part of the FL densitj*^. 
The cr-term does not contribute to the renormalization of 
since the cr-term can be eliminated in every step of 
the RG procedure by applying the transformation that 
does not affect the correlator Ro, as, discussed at the be- 
ginning of this section. That is why Y.^ and Ej, will be 
generated also for d < 4. Villain and Fernandez^ found 
from a non-perturbative RG that for d < 4 the defect- 
induced disorder flows under the RG to strong coupling. 
However, our study of the strong coupling limit in Sec- 
tion IV shows that this limit gives qualitatively the same 
result as the case investigated in this section. 

To summarize, the planar glass phase is characterized 
by (i) diverging shear and tilt moduli but a flnite com- 
pressibility, (ii) a transverse Meissner effect as well as 
a resistance against shear deformation, (iii) sample to 
sample fluctuations of the longitudinal magnetic suscep- 
tibility and (iv) an exponential decay of positional cor- 
relations in the direction perpendicular to the defects in 
d = 3. 

Since point disorder may formally become relevant be- 
low d = 4, we consider the stability of the planar glass 
with respect to weak point impurities. We find that the 
pinning energy due to point impurities in d = 3 behaves 



{np)^PoY^ / drV^p(r)e™'^°"e-("<^°)'<"'>/2 

- PQ^n,rapvlLyL,Lxe-^-'^^^-^ (63) 

and from this we conclude that weak point impurities are 
an irrelevant perturbation. Similarly, it can be shown 
that pinning energy of randomly distributed columnar 
defects decays exponentially with Lx and hence does not 
destroy the planar glass. 



D. Stability of the Bragg glass and the weakly 
pinned Bose glass 

In this subsection we continue discussion of the com- 
petition between pinning effects due to point impurities 
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and columnar and planar defects. We shall show that 
the weakly pinned Bose glass is stable with respect to 
weak point impurities but unstable with respect to weak 
planar defects. Moreover, we shall demonstrate that the 
Bragg glass phase is unstable with respect to both weak 
planar and weak columnar defects. The resulting phase 
diagram is shown schematically in Fig. [T] 

First wc discuss the stability of the Bragg glass phase in 
analogy to the test for stability of the planar glass in the 
previous subsection. We note that the correlation func- 
tions in the Bragg glass phase can also be obtained from 
a model with a uniaxial displacement field of FLsl. A 
uniaxial displacement field describes also properly charge 
density waves, a stack of membranes under tension and 
domain walls in magnets. Justified by these observations, 
we first consider a uniaxial displacement field in the di- 
rection perpendicular to the defect planes. At the Bragg 
glass fixed point in d = 3 the pinning energy of the planar 
defects behaves as 

^ p^^v-^U>lioL-^"l^\ (64) 

When we compare this energy to the pining energy of 
point impurities, 

Wp) ^ ^llJB*Afi) ^ L , (65) 

we find that planar defects are a relevant perturbation. 
Note that at the Bragg glass fixed point the system is 
isotropic, i.e., = Ly = = L. In Eq. (|65p we used 
the fact that the fixed point correlator on the length scale 
L behaves as R*p{0) - L-^ [cf. Eq. Similarly, the 

pinning energy of columnar defects at the Bragg glass 
fixed point scales as 

^^poV^^F^i^"'/'', (66) 

and drives the system away from the Bragg glass fixed 
point. 

A functional RG analysis of weak columnar defects in 
d = 5 — e yields a stable phase with a zero temperature 
fixed point that is characterized by a power law decay 
of the positional correlation function with an exponent 
Vc = ('''/3)^(5 — d) and a transverse Meissner effect^. 
One can expect that this phase, found for an uniaxial dis- 
placement field, applies to the case where the FL density 
is larger than the columnar defect density, corresponding 
to the so called weakly pinned Bose glass. In order to 
study the stability of this phase in d = 3 with respect 
to planar defects and point impurities, we compare the 
scaling 

^ PoVv'L^Ll/iDL-(^/'^\ (67) 
W^^^R*{Q)L'Ll^L,, (68) 

WP) - Po^^l^w^i-^"/'^' (69) 



of the different pinning energies, where we used L^: = 
Ly = L and columnar disorder fixed point correlator 
i?^(0) ^ L~^. We conclude that weak point impurities 
are irrelevant but weak planar defects are relevant in the 
weakly pinned Bose glass. 

Next we examine the stability of the Bragg glass phase 
by considering a vector displacement field. This displace- 
ment field reveals the triangular lattice structure of the 
FLL. By changing the orientation of the defect planes, 
the number of FLs that are pinned by the defects changes. 
In d = 3 wc have 

(77^ ^ po^v^Lyeoe-'^o'm/^ ^ Li-s, (70) 
(Tfc) ~ pov/^i^e-«^^/2 ~ L'~^/^, (71) 
^L^R*p{0)-L, (72) 

where g is given by Eq. ([25]), |Go| = in/iV^a) is the 
shortest vector of the reciprocal lattice and t] is the ex- 
ponent of the positional correlation function in the Bragg 
glass phase. Weak columnar defects arc always a relevant 
perturbation, while weak planar defects are relevant only 
if they satisfy g < gc — 3/2, i.e., only if they arc paral- 
lel to the main crystallographic planes of the FLL. Here 
we neglected the influence of weak planar defects on the 
elasticity of the FLL. In fact, the defects lead to an addi- 
tional anisotropy in the elastic energy which is associated 
with a larger energy for deformations with nonzero dyU^ 
and dzUx- Through a renormalization of the elastic con- 
stants g is renormalizcd downwards. Therefore, stronger 
planar defects lead to an increased gc > 3/2, rendering 
additional orientations of defects relevant. However, it is 
likely that in the case of a finite density of parallel de- 
fect planes, the FLL will rotate to a position in which 
it reaches maximum overlap with the defects. Then the 
planar defects will be parallel to the main lattice planes 
and the Bragg glass is unstable. 



V. FINITE DENSITY OF STRONG DEFECTS 

In this section we consider the FLL with a finite density 
of parallel, randomly distributed defect planes that are 
aligned to the magnetic field. The defects are assumed 
to be sufficiently strong so that the Larkin length [see 
Eq. ([54]) ] is of the order of the mean defect spacing or 
smaller, < £d. In this case the defect potential can 
not be treated perturbativcly with respect to the elastic 
energy, and a new approach is required. Here we derive 
an effective Hamiltonian that is defined only at the defect 
planes. We determine the ground state configuration of 
the FLL and calculate the positional correlation function. 
We show that a transverse Meissner effect as well as a 
resistance against shear deformations appear also in this 
case. 
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The Hamiltonian in rf = 3 dimensions reads 



= Ti-o 



(73) 



where Tio is the elastic Hamiltonian given by Eq. ([2]) and 
TioA is the pinning energy of the defect plane at the 
position X = Xi, see Eq. ([22]) . No denotes the number 
of defects. In this section we consider a simplified model 
involving only uniaxial displacements perpendicular to 
the defects u = ux. In Sec. I VIII we shall discuss the 
implications of the generalization to a two-dimensional 
vector displacement. The part oi Hd that describes the 
coupling of the pinning potential to the slowly varying 
part of the PL density (~ dxu) leads only to an increase 
in the PL density at the defects and can be eliminated 
by applying the transformation 



Mr 



u{r) - vpo/cii^Q{x ■ 

i=l 



(74) 



where Q{x) — 0, x < and Q{x) 
Hamiltonian then becomes 

Nn „ lt/i]G 

H = 7io-2wpo^ / ^ cos {kGD[u{x,y, 



1, a; > 0. The 



z) - a,:]}, 
(75) 



where = Xi +vpo{i — l)/cii. For simplicity we assume 
that all defects have the same strength. 

In order to obtain a Hamiltonian that is isotropic 
in the j/z-plane, we introduce the rescalcd coordinate 
z' = z ^ C66 / C44 and define u'iy^z') = u{y,z). We shall 
omit the primes below. We proceed by studying the 
ground state of the displacement field for a given distri- 
bution of planar defects, assuming that a strong defect 
potential suppresses thermal fluctuations. First we solve 
the saddle point equation in the gap between the defects 
with prescribed, but arbitrary, boundary conditions at 
the defects Ui(y,z) = u{xi,y,z) with Fourier transform 
Ui{q) = u{xi,q), H = {qy, Qz)- In the ground state con- 
figuration the FLs are completely aligned to the defect 
planes and u{xi, y, z) is independent of y, z. However, we 
derive the saddle point solution and the effective Hamil- 
tonian for a more general displacement field configuration 
at the defect planes since this will be necessary for a dis- 
cussion of the transverse Meissner effect below as well 
as the FL dynamics in Sec. IVII In the following we use 
the notation \Ai = Ai^i — Ai for any quantity A. The 
solution of the saddle point equation between two defect 
planes reads, with x €E 

■"(^' A ^ ^ smh [q (Xi+i - x)\ 



sinh (q'Axi) 
Ui+ijq) 
sinh {q' Axi) 



sinh[g'(a;-a;,)], (76) 



where q' = \J c^^jcwq. Note that we have solved the 
saddle point equation beetween the defects within a con- 
tinuum model and not on the lattice. This amounts to 
setting the momentum cutoff A — > cxd. After substitut- 
ing this solution into the Hamiltonian of Eq. ((75)) and 
integrating over x, the Hamiltonian reduces to 



U = 



VC11C44 ( d^q 



(2^)2 



Wd-1 



|uj+i(q) - Mj(q)|^ 




sinh ^gA 
+ (|M»(q)P + |u,+i(q)p) tanh ' 

-2upoW — ^ / ^ cosffcCn [ui(y, z) - aj]}. 

V C66 ^.^^ Jy,z 

(77) 

A similar Hamiltonian has been obtained for a Luttingcr 
liquid with point impurities^. 

Next we study the ground state of the FLL. In the limit 
w — > 00 the FLs are completely aligned to the defects and 
Mi(y, z) = Ui. Then the Hamiltonian of Eq. (j77p becomes 



9 2^ 



. , Axi 

Nd l^/i]G 
i=l k>0 



k— [Ui - at 



(78) 



where L is the system size. For Ui — irii + ai, where is 
an integer number, the energy gain of the FLs from the 
defect potential is maximal. We determine rii such that 
the clastic energy is minimal and find the ground state 
configuration to be degenerate and given by^ 



Ac 



(79) 



where n is an integer number. Note that this is 
the ground state configuration for an arbitrary defect 
strength in the special case when Aai/i — [Aa^/^Jg = 
for all i. Then the FLs are just shifted in order to gain 
energy from the defects without any elastic energy loss. 
However, for randomly distributed defect planes that sat- 
isfy £_D » £, Aai/£— [Aai/£]Q is uniformly distributed in 
the interval [—1/2, 1/2]. Using the central limit theorem, 
we find that the positional correlation function decays 
exponentially fast in the x direction. 



x-» £1 



(80) 



with ~ Qio/T^^- This shows that the limits of weak 
and strong planar defects lead to the same behavior of the 
positional correlations in d = 3. The correlation length 
in both cases is determined by the Larkin length of the 
defects. 
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= (n + 1)1 



FIG. 4: Schematic illustration of a droplet with radius R and 
width w a,t a. defect plane. 

From the shifted boundary conditions z — > oo) = 
M° + £ and Uiiy, 0) = one can obtain also the interface 
tension and it turns out to be finite. We do not quote 
the result here since it is cutoff dependent and hence non- 
universal. Also, a general expression that is valid for all 
ratios of the elastic constants is not available. A simi- 
lar analysis shows that the surface tension Ej, is finite. 
Hence, strong defects lead to a transverse Meissner effect 
as well as a resistance against shear deformations. The 
finite values of and T,y might be interesting to probe 
exp eriment ally. 



VI. FLUX LINE CREEP 

The technologically the most interesting property of 
type-II superconductors is their ability to carry a bulk 
current with as little dissipation as possible. The Lorentz 
force acts on FLs and hence gives rise to dissipation^. 
Pining centers play an important role in preventing FL 
motion and lead to a nonlinear resistivity Inp ^ — 
for J ^ which depends on the so-called creep expo- 
nent fj^i^. In this section we study the effect of planar 
defects on the FL dynamics in the direction perpendicu- 
lar to the defect planes. The defects are assumed to be 
parallel to the applied magnetic field. For a single defect 
plane we show that the creep exponent is /i = 1, apart 
from logarithmic corrections. We find that many pla- 
nar defects act as a more effective source of pinning than 
point impurities^ and columnar defects^ i^^i^^ . They con- 
siderably slow down the FLs in comparison to the Bragg 
glass and the Bose glass, leading to a creep exponent 
/i = 3/2 for the planar glass. 

A. Single defect 

First we discuss the creep of FLs in the presence of 
a single planar defect, aligned to the applied magnetic 
field and without point impurities, for currents parallel to 
the defect and perpendicular to the magnetic field. The 
Hamiltonian is given hy TL — Ti-o +Ti.D + force, where 



the part describing elastic deformations, is given by 
Eq. the defect pinning energy Hd by Eq. ([22]) and 
the Lorentz force contribution reads 

H/orce = - y" d^r {J(r) xB(r)}-u(r). (81) 

J(r) is the current density, B(r) is the magnetic induc- 
tion and the speed of light is set to one. Since we will 
consider only FL dynamics normal to the defect and since 
the defect potential depends only on the perpendicular 
displacement field, it appears plausible to simplify our 
model and to consider a uniaxial displacement field in 
the direction perpendicular to the defect plane. 

The defect plane is a relevant perturbation in the Shub- 
nikov (mixed) phase for all orientations of the plane (in 
contrast to the case when point impurities are present) 
since thermal fluctuations do not roughen the FLL in 
d = 3. This can be easily seen from the RG equation 
p8|) by setting g = 0. Under renormalization the weak 
defect potential that couples to the periodic part of the 
FL density grows and flows to strong coupling. Since we 
are interested in the small current densities J ^ which 
probe large length scales, we shall study only the strong 
defect plane below. Our results apply also to weak de- 
fects at sufficiently large length scales or small currents. 

In the absence of a current, the FLs arc aligned to 
the defect plane and the ground state is highly degen- 
erate. Different ground states differ by a shift of the 
displacement field by i because the energy does not de- 
pend on which FLL plane is pinned by the defect. This 
degeneracy is broken when a current is turned on and 
the original ground state becomes unstable. The system 
now evolves into a new metastable state that is lower 
in energy and in which the FLs are shifted by £. This 
process is enabled via the formation of droplets which 
are nuclei of new metastable states. The competition 
between bulk energy gain and elastic energy loss deter- 
mines the energy (and size) of the critical droplet. The 
energy of the critical droplet corresponds to the energy 
barrier that the FLs have to overcome when evolving to a 
new state. Thermally activated FL hopping over barriers 
with energy E^^^p(J) determines the resistivity through 
the Arrhenius la w^^i^^ 

p{J) ^ e-^--('^)/^. (82) 

Therefore, we need to estimate E*^^^j^{J). 

We proceed by deriving an effective Hamiltonian that 
is defined on the defect plane (x = 0). By integrating out 
the displacement field off the defect we get (see Appendix 

El) 

We// = i j ^|u(0,g„g,)ng(0,q)]-i 

„ G 

-2vpQ / dydz ^ cos [fcG_Dw(0, y, z)] 
fe>o 

- u(0, c[ = 0) [ dxJ{x)B{x) . (83) 
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Here q = {qy,qz) is the in-planc momentum, 
m(0, qy, q^) = u{x = 0, qy,qz) and 



arctan 



e(0,q') = 



7i"|q'|VciiC66 



(84) 



where q' = (qy, \fc^Alc^qz) and A = In jX. Since the 
system is translationally invariant in the yz-plane, the 
current density and magnetic induction depend only on 
X. In order to simphfy the computations, we make the 
system isotropic in the yz-plane by the rescahng z' = 

1/2 

•2(c66/c44) and ■u'(r') = u[v). In the following we will 
omit the primes. 

The critical droplet is a solution of the saddle point 
equation for the Hamiltonian of Eq. ((83| with fixed 
boundary conditions u(0, p ^ 0) = (ri + 1)£ and u(0, p 
oo) = nl with p = [y^ + z^)^/^ and n integer. For the 
precise solution see, e.g., a related discussion for a sin- 
gle strong impurity in a Luttinger liquid by Giamarchi— . 
The shape of the droplet is characterized by its radius R 
and the width w of the droplet wall so that the displace- 
ment field obeys approximately 



ii(0,p) 



(n+1)^, p€ [0,i?] 

nt, p e [i? -t- w, oo]. 



(85) 



The exact shape of the droplet wall is not essential for 
the discussion that follows. We assume that it smoothly 
interpolates between (n -I- l)i and ni. The width w of 
the droplet wall does not depend on the radius for small 
currents J ^ 0. In this limit the critical droplet radius is 
much larger than the width w so that the energy loss 
is balanced by the energy gain from the Lorentz force (see 
e.g.— and references therein). The critical droplet radius 
and energy is determined by maximizing the droplet en- 
ergy Edrop{R)- 

Since across the droplet wall the FLs are not aligned to 
the defect plane, a strong defect tends to reduce the width 
of the wall. Then, the large q behavior of the propagator 
ofEq. ([84]) . [^(0,q)]^^ « (7rc66Q^)/A, becomes important 
at the wall since it describes elastic deformations on small 
length scales. For sufficiently strong defects [w ~ A~^) 
the precise form of the droplet wall is determined by the 
interplay between the elastic energy ~ q^ and the defect 
pinning energy. The energy loss for FLs at the defect 
plane is^ Ecore ~ Rv^^^, as known from droplets in the 
sine- Gordon model. 

The elastic energy loss outside the plane, due to the 
deformation of the FLs at the defect plane is determined 
by the low q behavior of the propagator of Eq. ([84|) . 
[5(0, q)]^^ ~ 2^/ciiC66|q|, since the deformation oc- 
curs across the large scale R. It captures the three- 
dimensional nature of the FLL by its non-local form 
~ |q|. This is obvious when the elastic energy is written 
as 



'Hel 



= / d^n / d^r.M^ll^l^ 

47r 7 J (ri - Y2Y 



(86) 



Here ri and V2 lay in the defect plane and satisfy 
|ri — > A. The long ranged elasticity in the effec- 
tive two-dimensional elastic Hamiltonian of Eq. (j86p re- 
sults from fluctuations outside the defect plane that have 
been integrated out. Since R ^ w, the precise form of 
the droplet wall is not important for estimating Tipj and 
we can assume w = (corresponding to v s- 00). Then 
we obtain for the elastic energy 



Eel « 2V^II^ri?log (R/X) . 



(87) 



This result can be interpreted as the energy of charges 
of equal sign (corresponding to kinks in the displacement 
field) which are placed along a circle of radius R and 
interact via the three-dimensional Coulomb potential. 

The energy losses mentioned above are balanced by an 
energy gain due to the Lorentz force that is described by 
the last term of Eq. 



E 



force '■ 



iR^TTj— / dxJ{x)B{x) 

C66 J 



When we estimate E force, a finite width of the droplet 
wall can be neglected since R ^ w. The total droplet 
energy then reads 



Edr 



op 



Eel + Ec 



aR log 



E 



force 



R(3 - 7i?2 



(89) 



where a = 2(011044)^^^^^, /? ~ v^/'^ and 7 = 
7r(c44/c66)^^^^ J dxJ{x)B{x). The creep rate is deter- 
mined by the droplet with the largest total energy Edrop 
which is called the critical droplet. Solving the equations 
duEdrop = and d\Edrop < we find the size R* and 
the energy E*^^^^ of the critical droplet. Increasing of the 
droplet radius beyond R* does not cost any energy and 
droplet freely expands. For R 00 the system reaches 
a new metastable state in which all FLs are shifted by £ 
perpendicular to the defect. The nonlinear resistivity is 
given by the Arrhenius law. In the limit of a vanishing 
current density, i.e., for large /3/7, we get 



P' 



exp ■ 



47T 



/3 + a log 



2A7 



(90) 



Prefactors are not determined here since in the limit of a 
vanishing current density the current-voltage character- 
istic is dominated by the exponential factor of Eq. (|90|1 . 
The result for p yields the creep exponent p — \ plus 
logarithmic corrections. To estimate the coefficient 7 we 
need to know how the current density and the magnetic 
induction vary in space. This is a tedious analysis which 
goes beyond the scope of the present study and is left for 
further investigation. 

Next we examine how randomly distributed weak point 
impurities around the defect plane affect the FL creep for 



17 



J 0. A weak defect is relevant in the RG sense only 
for 5 < 1 and it then flows to the strong coupling limit. 
Criteria for the relevance of a strong defect with an ar- 
bitrary orientation are not available. Therefore we study 
below only the strong coupling limit for a defect that is 
oriented parallel to the main FLL planes. We expect that 
the liberation from the defect plane is the limiting factor 
for the FL motion so that the creep exponent is reduced 
compared to its value in the BG phase. 

The shape of the droplet is again given by Eq. ([55]) 
but now the impurities control the fluctuations of the 
FLs outside the defect. Without point impurities, the 
displacement field decays outside the defect plane as 
u{x,p = 0) w {t/2){R/x)'^ for a; > i?. Point-Hkc impuri- 
ties induce additional displacement fluctuations and the 
droplet-induced deformations are no longer long-ranged. 
This can be seen from the correlation function in the 
presence of a relevant defect that follows from Eq. (|40|) , 



droplet is then of the form 



'u{x,p)-u{Q,p)f 



BG 



(91) 



The subscript BG means that the correlation function is 
computed at the BG fixed point. Since the right hand 
side of Eq. ((9T|) is of the order a? for x w La we con- 
clude that the droplet extends up to La from the defect. 
This yields for the energy gain from the Lorcntz force the 
rough estimate E force ~ JBLaR^i where J is the mean 
current density inside the droplet. Hence, in the limit 
J — > resistivity is 



p{J) ~ exp 



(92) 



where C'l and C'2 depend on T, B, on the strength 
and concentration of the impurities, and on the defect 
strength. This result shows that a single relevant defect 
plane indeed slows down the FL creep in comparison to 
the BG phase. 



B. Finite density of weak defects 

Here we consider FL creep perpendicular to many weak 
defects with random positions but in the absence of point 
impurities, cf. Sec. IIVI The motion of FL bundles un- 
der the influence of the Lorentz force is again driven by 
the nuclcation of critical dropletSi^. A typical droplet is 
schematically shown in Fig. \5\ For small currents, the 
droplet extends over many defect planes in order to bal- 
ance elastic and pinning energy loss with bulk energy 
gain from the Lorentz force. For small current densities, 
the FLL is properly described in terms of the interface 
tensions Ej, and S^, see Eq. (|5ip . which are appropri- 
ate on sufficiently large length scales. The energy of the 



E, 



drop 
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(93) 



The elastic energy cost for the formation of the droplet 
consists of two terms. The first term of Eq. ((93)) is the en- 
ergy of a wide domain wall of width ^ L^ parallel to yz- 
planc. The second term of Eq. ([M)) describes the energy 
of a narrow cylindrically shaped domain wall perpendic- 
ular to yz-plane. In the estimate of Edrop we have taken 
into account that the elastic energy and the energy from 
planar disorder scale in the same way. The last term 
of Eq. ((93| is the energy gain from the Lorentz force. 
J (B) is to be understood as the mean current density 
(magnetic induction) averaged over the defect spacing. 

Note that we have used again the rescaling z' = 
(c66/c44)^/^z. In Eq. ([M)) we have taken into account 
the logarithmic roughness of the displacement field, cor- 
responding to the roughness exponent C = 0- We have 
ignored logarithmic corrections. The a-term of Eq. (|45p . 
that is responsible for the roughness exponent ( = 1/2, 
can be eliminated by the simple transformation that was 
discussed in Sec. IIVI and hence it does not affect the FL 
dynamics. In case of a potential breakdown of the e ex- 
pansion (cf. Sec. lIVBp in d = 3 and the existence of a 
strong coupUng fixed point^ see Sec. IVI CI 

To determine the critical droplet we solve dL^Edrop = 
and djiEdrop = 0- We find for the critical radius R* 
and the critical length L* of the droplet 



R* 

S7/ 



l: 
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ciiP jBe 



(94) 



This yields for d = 3 the nonlinear resistivity in the limit 

J < Jd, 



piJ) 



Jo^C ^^^-'^''\ (95) 



Here C is positive numerical constant of order unity. Thus 
the non-linear resistivity is considerably reduced com- 
pared to the Bragg glass phase and to a single defect 
plane in the presence of impurities. 

In the similar way one can consider the creep in the 
presence of randomly distributed columnar defects that 
are aligned to the applied magnetic field and have a mean 
spacing that is larger than FL spacing. Based on a func- 
tional RG in d = 5 — e dimensions^, wc obtain the creep 
exponent /.( = 1. This result is in agreement with the one 
derived in Blatter et al.— by other means. The result is 
expected to apply to the weakly pinned Bose glass phase. 



C. Finite density of strong defects 

Here we discuss the analog of the previous subsection 
in the limit of strong defects, sec Sec. |Vl The ground 
state degeneracy given by Eq. ([75]) is broken when a 
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current is applied. The system evolves between differ- 
ent ground states via the formation of critical droplets. 
For small currents J ^ cul/ (Bi'jj) the critical droplet 
extends over many defect planes and we obtain the 
creep exponent fi = 3/2. For moderate currents with 
cii£/{Bel)) < J < w/(0oC^n) the droplet forms only 
at a single defect plane and we recover Eq. (j90|) with 

1/2 

7 w (C44/C66) UdJB, i.e., a creep exponent ^~\. 

We assume that the saddle point solution for Ui with 
fixed boundary conditions Ui{p — > 0) = uf'^^^ and 

Ui{p 00) = Uj-"^ obeys Eq. ([55]) . At each defect plane, 
the radius Ri and the center of the droplet can be differ- 
ent. However, it is plausible to assume that in the saddle 
point configuration the droplet is centered at the same 
lateral position in each plane and all droplets have the 
same radius R since the droplet tends to maximize its 
volume while keeping the surface minimal. Specifically, 
we assume that the droplet is located between the sth 
and the (s -f- m)th defect plane. 

The width of the droplet wall for a sufficiently strong 
defect potential satisfies w <^ £d and w <C R- The pre- 
cise form of the wall is not important in finding the en- 
ergy of a domain wall parallel to the defects as well as 
the bulk energy gain. Hence, we set the width of the 
droplet to zero which yields for the Fourier transform of 
the displacement at the ith defect plane 



w,(q) = 27r Ri 



MqR) 



+ 4'\2nrS{<i), (96) 



where Ji is the Bessel function of the first kind. Due to 
the transport current we have to add to the Hamiltonian 
of Eq. ([77| the additional energy 



Nd-1 



H 



force 



+ 



Ax., 

Ax. 



(x-x,)}, (97) 



where /(x) = (c44/c66)^^^ J(x)_B(x). By substituting 
the Eq. dMl) into Eq. ^7}, we find that the droplet 



energy has a different form for R iDy/cee/cn and 
R<.iD ^cee/cii. 

First we discuss the case R ^ v^^ee/cii • The 
droplet energy reads then 

Edrov = Sx-R^TT + S„2i?7rm^D - JBiR'^T:m£D\ — ■ 

V C66 

(98) 

In order to explain and interpret the first term of 
Eq. ([M)) . we consider an excited state with 



(nil) r- 

u) for 



for i < k 
i > k . 



(99) 



This state describes a domain wall parallel to the defects. 
Using Eq. ([75)) . it can be shown that for a given disorder 



■2R 




FIG. 5: Scliematic representation of a cylindrically shaped 
droplet of radius R and lengtli Lx that extends across more 
than one defect plane. It drives the system into the new 
metastable state -I- £. Note that domain walls parallel to 
the defects are wide (not shown here), while the cylindrically 
shaped domain wall is narrow (of width w) . 



realization the energy cost of such wall per unit surface 
area is 



S]J(fc) = cii 




1±2 



Aak 



(100) 



Since the droplet consists of two such walls (sec Fig. [5]), 
the surface tension in the first term of Eq. ([55]) is given 
by 



- v+ 



Sj(s-l) + E-(,s + ,n), 



(101) 



The droplet takes advantage of fluctuations in the sur- 
face tension 'E^ and the lowest value of "E^ for a 
drop let of le ngth = mio is typically of the order 
ci 1^/044/ C66^^ / Lx when £d ^ £■ We point out that this 
result is in agreement with the statistical tilt symmetry. 

The second term of Eq. ([98]) is the energy cost for the 
domain wall perpendicular to the defects with surface 
tension . We do not provide here an explicit expression 
for Ej,, for the reasons discussed in Sec. |Vl Note that 
carries information about the strength and density of 
defect planes. By ignoring the spatial variations of J and 
B, we get the average bulk energy gain to be given by 
the last term in Eq. ([M]) . 

By solving dmEdrop = and d^Edrop = we deter- 
mine the radius R* and length L* = m*£D of the critical 
droplet. With this parameters, Eq. yields the en- 
ergy of the critical droplet. We find that this energy 
leads, up to unimportant prcfactors, to the same nonlin- 
ear resistivity as in the weak pinning case of Eq. (j95p . 
This result is valid for sufficiently small currents such 
that R* » £D\/cee/cii and m* ^ 1. The latter con- 
dition results from treating m as a continuous variable 
which is a reasonable approximation for critical droplet 
that extend across a large number of defect planes. These 
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conditions translate to the requirement J <C cii^/(J5£|,) 
for the current density. 

Next, we shortly discuss the droplet expansion when it 
reaches the radius R* and the length L*. By analyzing 
the Hessian matrix of Edrop we find that the eigenvector 
that points into the direction of the free droplet expan- 
sion has a cc-component that is much smaller than its 
p-component, i.e., an expansion of the cross-section of 
the cylinder is favorable over an expansion of its length. 
In order to describe a potential growth in the x-direction, 
one has to know the set of numbers that depends 

on the disorder realization. However, the droplet will get 
stuck between planes with low T.^ and a further expan- 
sion along the x-axis costs energy. 

Droplets occur and expand independently across the 
entire sample. After the FLs have moved in some regions, 
the boundaries of these regions will be favorable sites for 
emergence of new droplets^. Indeed, after the formation 
and expansion of a droplet up to the (i — l)th defect plane 
such that Uk = uj^"^^'' for k < i — \ and it^ = uj^"' for 
k > i, the new surface tension reads 



^xW ~ cii- / — — 
V C66 «n 

= -s-«<o. 













Jl 




. 27r 





(102) 



The reason for this result is that after the droplet expan- 
sion an additional FLL plane appears with respect to the 
initial ground state between ith and (i — l)th planes and 
the FLs are compressed. The formation of a new droplet 
at the zth plane is favorable because it allows the system 
to relax into the new ground state configuration between 
{i — l)th and ith plane. Since the droplet described by 
Eq. ([98|) has the longest life time, the resistivity is deter- 
mined by p ~ exp (y—E*i^^^/T^ , where the critical energy 

^drop given by Eq. evaluated at R* , L*. 

When comparing Eq. ([95|l with the creep exponent 
11= 1/2 of the defect free BG phase, we see that defect 
planes act as a more efficient source of pinning in stabi- 
lizing superconductivity than point impurities. However, 
we have considered only typical droplets in estimating 
Yjx- For system sizes L ^ £d it is likely that rare regions 
with untypically large T,x will appear and in turn deter- 
mine the resistivity. We leave this problem for further 
investigations. 

Next, we consider the case R ^ f^i^cgg/cn. From 
the second term of Eq. (|77p we find that for m defect 
planes the energy loss given by Eq. (jST]) . The first term in 
Eq. ([77| , that describes the coupling between neighboring 
defects, provides a much smaller contribution than the 
second term and can be neglected. Then the defects are 
effectively decoupled, and the nucleus energy is 



Edrop — 'n'iE, 



single ■ 



(103) 



Here Esingie is the energy of the droplet that appears in 
the case of a single defect plane. It is given by Eq. ([M]) 
with the system size replaced by the mean defect distance 




FIG. 6: (Color online) Array of edge dislocations that is lo- 
cated at the defects (thick blue lines) in order to relax shear 
strain. The Burgers vector is parallel to the defects. 



to- The nucleus energy grows with increasing m so that 
it is minimal for m = 1 and the critical droplet is located 
at one defect plane only. Nonliner resistivity is then again 
given by Eq. but with 7 « {ca/c&Qf^UDJB. This 
result is valid for intermediate currents cut/ [Bt^j] <C 
J<z;/(M^i,). 



VII. DISCUSSION 

For a finite density of randomly distributed parallel 
planar defects with the magnetic field aligned to them, 
and with a mean defect spacing that is larger than FL 
spacing, we find a new phase of FLs at low tempera- 
tures, the planar glass. We considered mainly a simpli- 
fied model with an uniaxial displacement field (which is 
also applicable to a wide class of other systems). Here 
we comment on possible consequences of this simplifica- 
tion by taking into account also the displacement field Uy 
parallel to the defects. The part of the defect Hamilto- 
nian that describes the coupling of the defect potential to 
the slowly varying part of FL density can be eliminated 
by transforming described by Eq. ([74]) . For strong 

planar defects each defect plane is occupied by a single 
FL layer and hence ^^"'^(cci, y, z) = trii -\- ai for all y, z 
in order to maximize the pinning energy gain. Even in 
the absence of point disorder the displacement Uy does 
not vanish. This can be seen most easily in the case of 
isotropic elasticity where the following relations hold 



(jdxUx = -dyUy, 



Cll ~ C66 
Cll -I- C66 



(104) 



Here a is the Poisson number with — 1 < a < li^ The 
strain dxU^ in the gap between the defects at Xi+i and 
Xi is 

dxUx^l + ^+ip, (105) 

CiiAx; AXi 

where we used the notation AAi = Ai^i — Ai for a vari- 
able A. The difference of the strain dyUy in neighboring 
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gaps is then 
AdyUy « —a£ 



Axi+i 



Aril 
Axi 



1 



vpo 



Cut yAxi+i 



1 

(106) 



which is of the order -izal/iD- On the scale Ly this im- 
phes Auy ^ ±a£Ly/£u. To avoid a diverging shear en- 
ergy, dislocations with a Burgers vector parallel to the 
y-direction occur at the defect planes (see Fig.[S]). Their 
distance along the y-direction is of the order £d/o'- The 
energy of a pair of edge dislocations with anti-parallel 
Burgers vectors at a distance iu is^ 



77 ^ '^66 2 r , / 
J^edge ~ '2^'^ \~a 



(107) 



This energy has to be compared with the energy gain 
from the defects which is of the order 
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APPENDIX A: REPLICA HAMILTONIAN FOR 
THE DEFECT FREE SYSTEM 

In this Appendix the Hamiltonian for defect-free sys- 
tem is derived, using the replica approach for averaging 
over point impurities. The pinning energy of randomly 
distributed impurities reads (see Sec. |TT| 



E 



D 



-L^vioPo ■ 



(108) 



Hence for acQ^a^^ <C tnv the energy of the dislocations 
is overcompensated by the defect planes and dislocations 
will be present. 

In general, the network of additional FLL sheets 
spanned by the dislocations will be complicated. The 
network follows from the solution of the equations of 
two-dimensional elasticity with the boundary condition 
Ux{xi^y, z) = ulP''\xi,y, z) and the dislocation density 
by{xi,y) at each defect. The energy has to be minimized 
first with respect to by{xi,y) and the with respect to n^. 
The resulting state is completely ordered along the z- 
direction. It is also ordered in the sense that the interface 
tensions Ti^ and are non-zero. A change in the bound- 
ary conditions with Ux{x, y,z ^ oo) = Ux{x, y, z ^ 0)+i 
increases again the energy. Hence the transverse Meiss- 
ner effect as well as the resistance against FLL shear- 
ing perpendicular to the defects are still present. Bond- 
orientational order— persists since disclination pairs re- 
main bounded in the cores of the edge dislocations. 

Since the Burgers vectors of the dislocations arc al- 
ways parallel to the defect planes, creep along the x- 
direction is not facilitated. Under the assumption that 
the distribution of E^, is uniform even in the presence of 
dislocations we recover the creep law of Eq. ([55]) . To de- 
scribe creep parallel to the defects one has to take into 
account the interaction between dislocation, a problem 
not considered so far— i^. We leave this for further stud- 
ies. For weak pinning qualitatively the same behavior 
can be expected on scales ^ Ld, in particular if the 
defect potential flows under the RG to strong coupling. 
If the sample exhibits two orthogonal families of (non- 
intersecting) defects, long range order in the x and y 
direction is destroyed even without point impurities on 
scales larger than Lo- The creep is then limited by the 
slowest mechanism and hence Eq. ((551) likely to be valid 
for all current directions in the xy-plane. 



rf^r yp(r)po < -Vxu(r) 



•iG[x-u(r)] 



(Al) 



where x — {x,y). If the system is characterized by a 
roughness exponent displacements vary with the scale 
L as w L'', ( < 1. The elastic energy scales as 
^d-2-i-2C^ and the first and second term of Eq. (|A1[) scale 
as l(''~2+2C)/2 L'^/'^^ respectively. These simple scal- 
ing arguments show that the coupling of the divergence 
of the displacement field to the disorder potential is ir- 
relevant with respect to the clastic energy in d > 2 and 
the second term of Eq. (|Aip is relevant for d < 4. Since 
we are interested in the behavior on large length scales 
in d = 3 dimensions, we can neglect the first term of 
Eq. (|Aip . After performing the disorder average, the 
replicated pinning energy reads 



2 2 ^ 



2T 



10 n 

^ / <5ax-x')^(z-z') 

a 1 Jr.r' 



Q,/3=l 

G[x-u°(r)]+iG'[x'-u'^(r')] 



G.G'^O 

„,2„ „2 n 



2T 



E E 



-i/^n X^ 



+Xr-u''(x+Xr,z)] 



G't^O 



2 2 

VpriimpPo / „ix(G+G') 



E E 



2T 

g-i[Gu°(x,z)+G'u''(x,2)]x _^(-Q')^ 



(A2) 



Using the relative coordinate = x' — x and the fact 
that ^{(xj.) is nonzero only for jx^l < ^, we approximate 
the slowly varying displacement field as u'^ (x -I- x^ , z) « 
u''(x, z), where S^-i (G) is the Fourier transform of 5^(x). 
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Since the displacement field varies slowly on the scale 
of the FLL constant, the integral over x vanishes for all 
combinations of G and G' except for G = — G' when the 
oscillatory factor e*''(*^+*^ ^ becomes one. The replicated 
pinning Hamiltonian can now be written as 

1 " r 

= -2T ^ / i?p[u"(r)-u''(r)]. (A3) 

APPENDIX B: EFFECTIVE HAMILTONIAN ON 
THE DEFECT PLANE 

In this appendix we present a functional integral ap- 
proach for the derivation of the effective Hamiltonian on 



the defect plane for the case of a FLL with point impu- 
rities and a single defect plane. Since the pinning energy 
of the planar defect involves only the displacement per- 
pendicular to the defect, u^, we integrate out Ux outside 
the defect and Uy across the entire sample. The partition 
function can be written as 

Z= fv^ivn) fvuir)e-^l[5[ux{rD)-^{rD)], 

(Bl) 

where H = Ho + Tip + Hd- The displacement at the 
defect is constrained to be Uxirn) = '■P^'^d) and this con- 
straint is implemented by (5-functions in the functional 
integral. The defect is aligned to the magnetic field and 
r^i is given by Eq. (|20|). After averaging over point im- 
purities we get 



•' '' a,Tn 

V[ip'^(rD)]e (B2) 

where Hq is given by Eq. ([T6|) and J P[(y3"(r£))] = n"=i / T^f^i^o)- Using the function integral representation 

n 5K(r,3) -^"(r,,)] = / P[A"(rz,)] ^° (B3) 

a,TD 

of the (5-function, we obtain for the effective replica Hamiltonian the equation 



e T = e 

up to a constant. Here (. . denotes the thermal average with respect to Ti.. For this average we obtain 

where (u^(q)uf (q'))-Hn = T(27r)''^^/(q)(5(q + q'). The effective Hamiltonian reads 

^"// = -EE2«fc^o / cos{fcGz,[<5-^"(rz,)]} + i5^— ^ / d'^-iqv'"(q) Q~.]M) i'^)^ (B6) 
where Q(q) = Gxx{x ~ 0,q) and here q is the in-plane momentum. 

APPENDIX C: DENSITY OSCILLATIONS IN THE PRESENCE OF AN IRRELEVANT DEFECT PLANE 

In this appendix we analyze the thermal and disorder average of the FL density around an irrelevant defect plane 
in the presence of point impurities by perturbation theory in v. For the local density variations we have 

(Jp[r,u(r)]) ^ lim / P[u"] 5p[r,Ui(r)] e-^[^S+^- ^"^O], (CI) 
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where 5p[r, u(r)] = p[r, u(r)] — po and ui(r) is the displacement field with replica index a = 1. To the zcroth order 
we get 



(<5p[r, u(r)]> = lim (<5p[r, ui(r)])„5 = 0, (C2) 

since (u^(r))'H^ = oo. To capture the physics correctly, we have to calculate mean FL density at least to first order 
in V (see the discussion in Sec. IIIID 



n 



{Sp[r, u(r)]) = - /? lim ^ {Sp[r, Ui(r)] Hz,(u"))h. 

= /3 hm {(5p[r,ui(r)] Hc(u2))«„" - (<5p[r, Ui(r)] ^^.^(ui))^^ }. (C3) 



M2(r) is the displacement field with replica index a = 2. First, we obtain the average of the long wavelength part of 
the FL density. It can be shown that 

(VxU^(r)cos[Gz,(5-Gou"(rz,)])^^ = Vx{Tsin(Gz3,5) (cos [Gi,u°(ri,)])„. g"'^ {vd ~ r)\GD) 

-Tcos{GdS) {sm[GDu''{rD)])H^g"-^{rD-r)\GD)}, (C4) 

where r = (x, z) and ^"'''(r) is the propagator given by Eq. pT)l . Since (sin [G£)U"(r£))])-Hn = (cos [G£)U"(r£))])-Hn = 
this contribution vanishes. 

A finite difference between the expressions {Sp{ui) H_d(u2))-hj and (Sp{ui) Hd(ui))-hj appearing in Eq. (jC3p can 
result from the thermal part of the propagator that is diagonal in replica indices, 



lim(<5pK(r)]7^^(u'^))^„ = lim{^ / ^ e'°-(e-'° ""W cos {Gz,[5 - (rz,)]}) } 

= l™o / ^ ^'""'^ + } , (C5) 

-i([G-u°(r)±GD-u''(r£,)]^) „ 

where I± = e "o and J^^ denotes the integration along the defect plane. Analyzing (/_), 

we conclude that it is nonzero only for G = — Gd (G ~ Gd) and 

J±«e— ^ ^ . (C6) 



This yields 



•^nl - 1 



^'"'■''»^" cos ICM^ - sn (^V F (,r^ ) , (C7) 



T \ 2ttc J \47rc|x — 5| / ' 

where F{x) = X^J^i ^ 2g+K-2 • ^'^'^ very small temperatures the main contribution is 



The result captures the large scale behavior since it is valid on scales larger than L > La- Here vi denotes the effective 
defect strength measured on the scale La, and \x — S\ is the distance to the defect plane. Additional contributions to 
Eq. (jCSp . coming from the higher harmonics in TYd, are less important since they are proportional to the coefficients 

Vk at scale L = La and their amplitudes decay as |a; — ^"""^ with integer k > 2. 
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APPENDIX D: DENSITY OSCILLATIONS IN THE PRESENCE OF A RELEVANT DEFECT PLANE 



In this appendix we study the displacement correlation functions and average FL density profile for a relevant defect 
plane in the presence of point impurities. As shown in the main text above, on sufficiently large length scales pinning 
effects can be taken into account through the boundary condition UxiYo) = at the defect plane. For simplicity we 
take the defect to be at the coordinate origin, i.e., we set (5 = 0. First we calculate the generating function 



Z«[j"(r)] 



T e 



E„/J°(r)u°(r) 



Using the representation of the delta-function of Eq. (jBSp we get 
where 



ri - r2 



7K 



(Dl) 



(D2) 



(D3) 



and Q is given by Eq. ((3T|) and Q"'^ is the inverse of given by Eq. ([T7)) . The displacement correlation function 
is given by the relation (u(r)u(r)) = lim„^o <5jtfr)fji(r) • Denoting by q the in-plane momentum, we get in 



momentum space 



j°=0 



jix,x;q) = lini 

n— i-0 



^^l' (0. q) - E (1^1 ' (q)^^S (1^1 ' -q) 



a7 



(D4) 



where Gpin.ij{x,x';r\\ — r'|) = T ^ {ui{r)uj{r')) with r = (x,r||). The indices i,j take the values x,y. All propagators 
and their inverse that appear in the previous equations have the form Xa.p ~ Sa.pX^ + Xn- The only nonzero 
contribution to the second term of Eq. (|D4p comes from the product of all "diagonal" parts (X^) of the propagators 
or only one " nondiagonal" (Xn) and two diagonal (in replica indices) in the limit n — s- 0. Denoting by Xij,a = 
Imin^Q {i\X a\j) , where a = d,n, one has 

Gpin,ij{x,x;q) = ^y,d(0, q) - ^„,d(|a;|,q) GjxAl^l^'i) Qd^i'i) 

+ ^jj,„(0,q) -^„,d(|a;|,q) §jxA\x\,<i) Qn\'i) + 2d^(q) GjxAl^l^'i) fe,d(|a;|,q) 
+ Qd'(q) g^xM,c^) gjx.M,(l)- (D5) 
For isotropic elasticity the relation 

g,xAHq) g^xAl^ln) Qd'(q) = (x-«)(x-j)^y,d(2|2^|,q) (D6) 



holds. After integrating the terms of the second and the third line of Eq. (jD5|) over q we find that the displacement 
correlations on scales larger than L„ ([38|) read 

(D7) 



(«i(r)u,(r)) = Ihn {Tgij(0,O) - T(x • • j)glli2\x\,0)} . 
Next, we shall calculate the disorder and thermal average of the FL density 



{6p{v)) ^poY. 



using 



and 



(e-»Gu) ^ Um(e'''*^"0 = lim e-3((G"iM)'> 
hm ([Gui(r)]2) = lim T [G' ^14(0) - (G • (2|x|, 0)] . 



(D8) 



(D9) 



(DIO) 
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Since lim„^o ^^^(0) is divergent we conclude that only terms with a reciprocal vector G perpendicular to the defect 
plane contribute in Eq. (jPSp and 



m>0 



{5p{r)) =2po}^ cos (mGux) { ] , (Dll) 



where m is an integer. 



APPENDIX E: SAMPLE-TO-SAMPLE 
FLUCTUATIONS OF THE MAGNETIC 
SUSCEPTIBILITY 

In this appendix we examine the influence of planar de- 
fects on the longitudinal magnetic susceptibility. An in- 
finitesimal change in the longitudinal magnetic field SH~z 
changes the Hamiltonian of Eq. (|42p in the case of an uni- 
axial displacement field perpendicular to the defects by 



SH 



4>aPo 
47r 



(El) 



Since the change of the magnetic induction is -B = 
PQff'odxU, the longitudinal magnetic susceptibility reads 



dSH, 



(E2) 



where F is the free energy. It is convenient to consider 
a generalization of this model to d dimensions where x 
is a d — 2-dimensional vector and xi is the component 
of X in the direction of the displacement u. Applying 
the transformation u —>■ u + hxi/cn, the additional term 
given by Eq. (jEip can be shifted away yielding 



-V + nD{u + hxi/cn) , (E3) 



nih,u) =no{u) 



2c 
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where F = L^Ly and h = SH^pQ(j)Q/{4'K). The pin- 
ning energy of planar defects Ti/^ can be written as 



Hoiu) ^ J drVDix)ps{u,r) + J drVD(x)pp(u, r) 

=nh+nl. (E4) 

where ps and pp are the slowly varying and periodic part 
of the FL density, respectively. Next, we would like to 



compute the average x- The free energy is given by 
F{h) = - Tlog Pue-^C*'")/^^ 

—V^Po—[drVDi^)-T\ogZi, (E5) 

311 Cll J 



2ci 



where Zi is the partition function for Tii (h, u) = Tip (w) + 
Hf^iu) + TPjj (u + hxi/cii). Using rephcas, the disorder 
averaged free energy can be written as 

h2 "yj^ _ Y 

F = V-T\im^ . (E6) 

2cii n^o n 

Here 'Zf = J 2?[u"] exp [-H^/T] where H'l is the replica 
Hamiltonian that follows from Ti.i(h,u). Since Ti.i{h,u) 
has the same statistical properties as 7i(0,M), i.e., it 
yields the same replica Hamiltonian, the only dependence 
on h in F comes from the first quadratic term in Eq. (jE6[) . 
Due to this so-called statistical tilt symmetry^, the dis- 
order averaged susceptibility 



X 



471 d^F {po^o) 



V dSm 47rcii 



(E7) 



is disorder independent. The important quantity are the 
sample-to-sample variations of the susceptibility. The 
free energy can be written as 



Fih) 



2ci: 



-V + Fo-T\og {e 



-Ho(u+hxi/cii)/T\ 



Ho J 

(E8) 



where Fq is the free energy of the system that is described 
by the elastic Hamiltonian only. To first order in pertur- 
bation theory with respect to ~ w we get 



AF(/i) = F{h) - F{h) = {Hd {u + hxilcii))^^ . (E9) 
For a system of linear size L^. in the x direction we find 



lf/^]G 

AF(/ii)AF(/i2) = 2{vpo)''npdiLyL,f / dx ^ e-("^-)'<"'>«o 

n>0 

I 



uGd 



(hi - h2) xi 



Cll 



(ElO) 



where Upd denotes the density of defects. Differentiation 
with respect to ft-i and h2 leads to 



where we have taken into account the irrelevance of ther- 
mal fluctuation. Since e > for d < 6, the sample-to- 



R'ii'mi (Lx 



L, 



(Ell) 
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sample fluctuations grow with the scale Lx- One cannot 
expect that this result is quantitatively correct for large 
Lx- However, qualitatively it demonstrates the relevance 
of defects and it is a sig nature of a glassy phase. For 
Lx > Lo we expect that (Ax)^/x^ approaches a finite 
universal value for d < 6. 



APPENDIX F: POSITIONAL CORRELATION 
FUNCTION FOR 6 > d > 4 

In this appendix the positional correlation function of 
the FLL with planar defects will be calculated, using 
perturbation theory and results from the functional RG 
analysis presented in Sec. IIVI The positional correlation 
functions have been calculated before for the FLL with 
point impurities for an uniaxial displacement field^ii and 
for a vector displacement field^i^. We perform the com- 
putations along the lines of these references. 

In a functional RG procedure, after integrating out fast 
modes in an infinitesimal shell with A/b < q < A, one 
can choose to keep the cutoff in momentum space fixed 
using the rescaling 



, X , z 

x = - z = — 
b 6x 



where z = {y,z). The displacement field is not rescaled 
due to the periodicity of Rd- This implies u{q) = 
5d-2+2xy'(-q') -y^c need to obtain the RG flow of the 
correlation function 



= Z-^ j Vu<{ci)u{cii)u{ci2) j Vu>{ci)e 



= z- 



X'it<(q)u(qi)M(q2)e" 



-;9H,(«<(q)) 



(F2) 



(Fl) 



where H = Ho + J Vd{x)p{u, r). Here u<(q) are modes 
that satisfy q < A/b and Tii is the Hamiltonian that ap- 
plies to the scale / = log& with b very close to unity. 
Using Eq. (jF2[) . we obtain a differential equation for 
(w(qi)w(q2)). We get 



i(qi)«(q2)) = fe2('^-™(7.'(q'i)u'(q'2)) (F3) 



where the only restriction on & is < A/fe, i = 1, 2. 



First, we calculate (u(qi)M(q2)) to lowest order in v, 
(u(qi)M(q2)) = (u(qi)u(q2))Ho + ^|im XI / ("i(qi)wi(q2)^£'K(x, zi) - u''(x, Z2)])^^ ^^^^^^ , (F4) 

where ui is the displacement field with replica index a = 1. We use the periodicity of Rd by writing R]j(u) — 
X)n '^^^ (uGdu). From this we find that at the planar glass fixed point 



-(27r) ^+^ 
(cii'zL)^ 



(u(qi)ii(q2)) = i R*D (0)<5(qi. + q2.)<5(qi.)5(q2.), (F5) 



where the rescaled fixed point correlator is i?^ (u) ~ [{u — i / 2)'^ ~ / 12\ for < u < i. Choosing b = 

A/ maxjgi^, (72x} in Eq- (|F3p so that it is justified to calculate the correlation function appearing on the right hand 
side of Eq. (jF3P at the fixed point, and using the result of Eq. (jFSp . we get 

(u(qi)u(q2)) = ^ efS{q,x + q2.)'5(qi.)5(q2.). (F6) 

Note that in Eq. ()F6P only displacements with qi^; = — q2a: are coupled as it would be the case for a quadratic 
Hamiltonian. Therefore one can write down an effective quadratic Hamiltonian in the q2;-momcntum space that 
reproduces the correlations to first order in e. The positional correlation function shows the power low behavior 

« lxr^(-/3)^ (F7) 

We point out that this result is valid only for d > 4. In d < 4 dimensions the part of the pinning potential related to 
the slowly varying part of the FL density also becomes relevant and further analysis is needed, see Sec. IIVI 
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APPENDIX G: LIST OF RECURRENT SYMBOLS 



V. 



'p 



Symbol Quantity Definition 

a flux line lattice constant 

B magnetic induction 

Cii i = 1,4,6 elastic constants Eq. ^ 

G reciprocal lattice vector 

g = |'?(f )^ parameter controlling the relevance of the single defect plane Eq. (|25|) 

Gd = ^ shortest reciprocal lattice vector perpendicular to the defect(s) Sec. IIIII 

TCq elastic energy of distortions of the ELL Eq. ([2]) 

Tip pinning energy of point impurities Eq. ^ 

Hq effective quadratic replica Hamiltonian for defect free system Eq. 

Hd pinning energy of planar defect Eq. (|22p 

J current density 

Larkin length Eq. pUjl 

La positional correlation length Sec. |IT] 

in mean distance between defects Sec. IIVI 

Ld Larkin length for planar defects Sec. IIVI 

riimp density of point impurities 

n]j unit vector perpendicular to the defect plane 

r/3 position vector of the defect plane Eq. (I20p 

Rp point disorder correlation function Eq. p4p 

Ru planar disorder correlation function Eq. (j46p 

R*j^ fixed point value of Rpi 

Sg (r) positional correlation function Sec. 

Sd surface of d dimensional unit sphere 

T temperature 

z) = u{xi, y, z) displacement field at the planar defect with x = Xi 

Ux displacement field perpendicular to the defects 

strength of point impurities Sec. ((IH) 



Vp pinning potential resulting from point impurities Eq. ([II| 

V defect strength Sec. |in 

Vd pinning potential resulting from planar defects Sec. IIV 

[x]g the closest integer to x 

6 defect distance to the origin 

C roughness coefficient Sec. [ll] 

r] power law exponent of positional correlation function in the Bragg glass regime Sec. [ill 

A momentum cutoff Sec. |IT] 

A penetration depth Sec. 

^PF roughness exponent in random force regime Sec. Hi] 

£,PM roughness exponent in random manifold regime Sec. |IT] 

£,BG roughness exponent in the Bragg glass regime Sec. 

^ superconductor coherence length 
correlation length 

p(u, r) flux line density Eq. © 

po background FL density Eq. jB]) 

^y(z) interface tension of domain wall parallel to x and z (y) axes Eq. ([5T 

cj)Q flux quantum Sec. [IT] 
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